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Abstract. Extending the approach of Grillakis-Shatah-Strauss, Bronski-Johnson-Kapitula, and 
others for Hamiltonian systems, we explore relations between the constrained variational problem 
mm X:C(x)=c c o € R r , and stability of solutions of a class of degenerate "quasi-gradient" sys- 

tems dX/dt — —M(X)X7£(X) admitting constraints, including Cahn-Hilliard equations, one- and 
multi-dimensional viscoelasticity, and coupled conservation law-reaction diffusion systems arising in 
chemotaxis and related settings. Using the relation between variational stability and the signature 
of dc/duj £ R rxr , where c(u) = C(X^) € R r denote the values of the imposed constraints and 
u> € W the associated Lagrange multipliers at a critical point X£, we obtain as in the Hamilton- 
ian case a general criterion for co-periodic stability of periodic waves, illuminating and extending 
a number of previous results obtained by direct Evans function techniques. More interestingly, 
comparing the form of the Jacobian arising in the co-periodic theory to Jacobians arising in the 
formal Whitham equations associated with modulation, we recover and substantially generalize a 
previously mysterious "modulational dichotomy" observed in special cases by Oh-Zumbrun and 
Howard, showing that co-periodic and sideband stability are incompatible. In particular, we both 
illuminate and extend to general viscosity /strain-gradient effects and multidimensional deforma- 
tions the result of Oh-Zumbrun of universal modulational instability of periodic solutions of the 
equations of viscoelasticity with strain-gradient effects, considered as functions on the whole line. 
Likewise, we generalize to multi-dimensions corresponding results of Howard on periodic solutions 
of Cahn-Hilliard equations. 



1. Introduction 

In this paper, we investigate from a general point of view stability of standing spatially peri- 
odic solutions of a broad class of equations with gradient or "quasi-gradient" structure including 
Cahn-Hilliard, viscoelasticity, and coupled conservative-reaction diffusion equations. Our main 
goals are to (i) understand modulational stability of periodic viscoelastic solutions with physically 
relevant viscosity/strain-gradient coefficients and multi-dimensional deformations, about which up 
to now very little has been known, and (ii) understand through gradient or related structure the 
commonality of results that have emerged through direct Evans-function techniques in this and 
other types of systems. 

In particular, we seek to better understand an apparently special modulational dichotomy that 
was established for a special model of one-dimensional viscoelasticity by Oh-Zumbrun in |OZlj . 
but later shown by Howard [Holj to hold also for the Cahn-Hilliard equation, asserting that, under 
rather general conditions, co-periodic and sideband stability are mutually incompatible: in particular, 
periodic solutions are always modulationally unstable. Both of these results were obtained by 
direct Evans-function computations, leaving the underlying mechanism unclear. Moreover, in the 
viscoelastic case, the technical condition of monotonicity of period with respect to amplitude left 
open the possibility that there might be special parameter regimes in which periodic waves could 
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be stable. Further, the method of proof does not appear to generalize to multi-dimensional planar 
periodic waves, as studied for example in [Yl IBYZj . 

Thus, these prior results left intriguing questions regarding the possibility of stable periodic vis- 
coelastic solutions. We note that, even though second-order reaction-diffusion equations cannot, 
by Sturm-Liouville considerations, support stable periodic waves, nor, for the same reasons, can 
conservative scalar Cahn-Hilliard equations, the Swift-Hohenberg equations give a well-known ex- 
ample of a higher-order reaction-diffusion equation that does support stable waves \Mi\ ISchnj , and 
generalized Kuramoto-Sivashinsky equations are another well-known example showing that con- 
servation laws can likewise support stable periodic waves; see jJNRZll IBJNRZ| IFSTj . Thus, the 
conservative structure of the equations of viscoelasticity is not a priori an obstruction to stability. 

A clue noted in |OZH IZ3| regarding possible common origins of the modulational dichotomies 
for Cahn-Hilliard and viscoelasticity equations are that both serve as models for phase-transitional 
elasticity, and both share a nonincreasing Cahn-Hilliard/van der Waals energy functional 8. Indeed, 
the Cahn-Hilliard equation is well-known to be the gradient system with respect to H^ 1 of the 
functional £, that is, 

(1.1) {d/dt)X{t) = -MV£{X(t)), 

where M = (id x )(id x )* = —d x is symmetric positive semidefinite. It has been conjectured |OZlllZ3] 
that the similar results obtained for viscoelasticity are due to this "hidden variational structure." 

Here, we show that this is true, in an extremely strong sense. Namely, we show that the equa- 
tions of viscoelasticity with strain-gradient effects may be put into the same form (jl.ip . but with M 
only positive semidefinite (i.e., M + M* > 0) and not symmetric- a structure that might be called 
"generalized gradient" or "quasi-gradient" form. More, we show that, under an additional property 
(|Eip of viscoelasticity, satisfied automatically in the gradient case M = M* , all of the results that 
have so far been shown for this system, including the results on modulational dichotomy, may be 
deduced entirely from this quasi-variational structure, constrained variational stability considera- 
tions like those used in [GSS] IBJKj and elsewhere to study stability of Hamiltonian waves, and a 
novel connection observed between the Jacobian and Lagrangian multipliers arising in constrained 
variational theory and the Whitham modulation equations formally governing sideband behavior. 
The latter we regard as the main new contribution of this paper. 

Indeed, we are able to put each of the mentioned applications into this framework and obtain 
a suite of similar results by these common methods. As in any investigation overlapping with 
classical calculus of variations and Sturm-Liouville theory, many of our results in the simpler cases 
considered here can be obtained by other means, and a great many have (no doubt more than the 
ones we have been able to cite) . Still more of our results can be obtained alternatively in a slightly 
weakened form by direct Evans function computations, which give a sort of mod two version of 
Sturm-Liouville applying in more general circumstances; see for example, the discussion of Sections 
HXT1 and 1231 

However, a significant number of our results appear also to be new, in particular (surprisingly) 
modulational instability of multiply periodic Cahn-Hilliard solutions (though see Remark 15 .71 on the 
two-dimensional case), Theorem 15.61 modulational instability of singly and multiply periodic cou- 
pled conservation-reaction diffusion (in particular, Keller-Segel) solutions, Theorem 15.91 universal 
modulational instability, without the technical condition of monotonicity of period vs. amplitude, 
of one-dimensional periodic viscoelasticity waves, Theorem 5.1; and modulational instability of the 
branch of planar multi-dimensional deformation solutions of viscoelasticity connected by transversal 
homotopy to the decoupled, one-dimensional case, Corollarv l5.2i By related techniques, we recover 
also versions of the results of |GZ} \Z3\ IPS2] on instability of homoclinic "spike" type solutions on 
the whole line; see Section 16.21 
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Also, though many of these ideas appear elsewhere in disparate sources, we feel that there is 
some value in putting a broad variety of problems into a common framework, continuing and 
extending the unification initiated in |PSll IPS2j . In particular, we both formalize the heuristic 
relation pointed out in |PS2j between coupled conservative-reaction diffusion and Cahn-Hilliard 
equations and extend this to the quasi-gradient viscoelastic case. We discuss also the natural 
connection to the Hamiltonian case of systems (d/di)X(t) = JV£(X(t)) with J = —J* skew 
symmetric. Here, interestingly, the appearance of conservation laws for generators —MA with M 
noninvertible is paralleled by the appearance for J noninvertible of conservation laws associated 
with noninvertibility of J rather than Noetherian symmetries; see (|3.7|) . Finally, in Section 15.51 we 
give a unified explanation of the previously-observed modulational dichotomies from the variational 
point of view, showing that the mechanisms lying behind these are unconstrained instability of A, 
together with a relation between Morse indices of the Floquet operators associated with MA and 
A, and the aforementioned property (lElj) that the center subspace of MA consist of ker MA. 

1.1. Discussion and open problems. As concerns our primary interest of classifying stability of 
periodic viscoelastic patterns, we have resolved some of the main problems pointed out in [OZl} lY]. 
but also pointed the way to some interesting further directions for investigation. In particular, it 
is an extremely interesting open question whether there could be a bifurcation to a stable branch 
of planar multideformational patterns as solutions are continued from the decoupled, effectively 
one-dimensional deformation, case. Another interesting open question is whether there are there 
other periodic viscoelastic solutions that are 'genuinely multi-dimensional" in the sense that they 
are not continuable from the decoupled case. These and related problems for homoclinic solutions 
are being investigated in [BYZj by numerical Evans function techniques. 

We note that an efficient alternative Evans method is proposed here for self-adjoint problems, 
based on the Sturm-Liouville duality between spatial and frequency variables x and A; see the 
last paragraph of Section 14.1.11 This also would be interesting to explore. Finally, it would be 
interesting to look for stable multiply periodic viscoelasticity solutions, perhaps corresponding to 
waves seen in buckling of viscoelastic materials [HLFY1 IBCS] ) ; see Section 15.1.11 

Note. We have been informed after the completion of this work that our analysis in Section 14.11 
overlaps with unpublished results of Bronski and Johnson on nonlocal Cahn-Hilliard equations. 

2. Preliminaries: variational stability and Sylvester's law of inertia 

We begin by recalling a fundamental relation ( (12.12|) below) underlying the Hamiltonian stability 
results of |GSSj . [B JKj . etc., between stability of constrained critical points and the Jacobian dc/doo 
of the map from Lagrange multipliers to constraints. This principle seems to go far back, at least 
implicitly, in the physics (and particularly the thermodynamics) literature [Cl ILL] . For a general 
treatment from the Calculus of Variations point of view, see for example [M] and references thereinQ 
Here, we give a particularly simple derivation^ based on Sylvester's law of inertia. 

2.1. Constrained minimization and Lagrange Multiplier formulation. Consider the gen- 
eral constrained minimization problem 

(2.1) min £(X), C = (Ci, . . . ,C r ) : H — > W, 

X:C(X)=c 



1 Curiously, the Jacobian (dc/duj)(u)*) emphasized in [GSS] and later works is never explicitly mentioned in [M] . 
even though the physical importance of the Lagrange multiplier is stressed. Rather, it appears in the equivalent form 
—y T A u} *y = — (VC, y), A^*y := VC, arising naturally in the analysis (and in practice most easily computable). 

2 Essentially a redaction of the argument of [M] 

3 



on a Hilbert space H. Here, the constraint function C is C 2 on H. Also, we assume that £ 
is differentiable on Diff(£) a dense subspace of H, and that for every X G Diff(£), V£(X) and 
V 2 £(X) are well-defined at least on a subset of Diff (£) that is dense in H. In addition, we assume 

(HI) VC is full rank in the vicinity under consideration. 

Under assumption (jHip . the feasible set {C(X) = Co} by the Implicit Function Theorem defines 
a local C 2 manifold of co-dimension r in the vicinity of any point X* with C(X*) = cq, tangent 
to VC(X*) ± J3 on which one may in principle carry out standard Calculus of Variations on local 
coordinate charts. 

Alternatively, one may solve an optimization problem on an enlarged flat space, following the 
method of Lagrange. Introduce the associated Lagrange multiplier problem 

(2.2) min S U (X) where £ U (X) := £{X) + u ■ (C{X) - c), 

w G M r are Lagrange multipliers, and the constrained Hessian (in X only) is given by: 

(2.3) A, := V 2 £u{X*) = V 2 £(X*) + u ■ V 2 C(X*). 
Critical points (X*,u*) G Diff(£) x W of §23$ satisfy C(X*) = c and 

T 

(2.4) V£(X*) = -^wjVC^r). 

3=1 

We recall the following well-known characterizations of stability. 

Proposition 2.1 (First derivative test). An element of X* G Diff(£) is a critical point of (|2.ip if 
and only if there exists uj* such that (X*,u>*) is a critical point of (12.2p for some lo* G W. 

Proof Critical points of flU) satisfy C(X*) = c and {V£(X*),Y) = for Y G VC(X*) ± , yielding 
V£(X*) G Span{V(7,(X*) : j = 1, . . . ,r}, or (J33D. □ 

Proposition 2.2 (Second derivative test). Assume hypothesis (IHlj) , Then the following assertions 
hold true: 

(i) A critical point X* G Diff(£) is strictly stable in Diff(£) (a strict local minimum of (|2.ip 
among competitors in Di&(£)) if A u * is positive definite on \/C{X*) 1 - nDiff(£), 

(ii) It is unstable in Diff(£), (not a nonstrict local minimizer of f)2. 1 1) among competitors in 
DiS(£)), hence in H, if A^* is not positive semi-definite on VC(X*)- L n Diff {£), 

Proof. Equivalents, (d 2 £ (X(t))/dt 2 )\ t=0 < (resp. g 0) for any path X(t) G Diff(£) with X(0) = 
X*, C(X(t)) = cq. Taylor expanding, we have that y(t) := X(t) - X* = tY x + t 2 Y 2 + 0(t 3 ), where 
Fi G Vclx*^ and {VC j {X*),Y 2 } = -^(Y 1 ,V 2 C J (X*)Y 1 ) for any j = l,...,r. Using, flM} we 
obtain that 

l -{d 2 £{X(t))/dt 2 )\ t= , = (Yt, lv 2 S{X*)Y x ) + (V£(X*),Y 2 ) 



{Y^-VH^Y,) -j^u* (VC 3 (X*),Y 2 ; 



1 



1 



= ^(y 1 ,^y 1 ), 

from which both assertions readily follow. □ 



We use this notation to denote {VC_, (X*) : j = 1, . . . , r} , slightly abusing the notation. 
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Remark 2.3. Note that the second variation conditions for the Lagrange multiplier problem are 
with respect to the constrained Hessian, and along the tangent feasible space VC(X*)- L only. 

Remark 2.4. The conclusions of this section are strongest, evidently, when Diff(£) = H, in which 
case strict stability implies variation in energy controls variation in H, or 

(2.5) ||Xi -X*|| 2 H < £(Xx) -£{X*) for ||Xi - A"*|| H sufficiently small. 

This holds, for example, for the Cahn-Hilliard and generalized KdV energies considered in Sections 
14.11 and 14.41 However, the Keller-Segel and viscoelastic energies considered in Sections 14.21 and 14.31 
are not even continuous on H, hence the more general phrasing here; see Remarks 14.131 and 14.121 

2.2. Sylvester's Law and relative signature. Let Q : dom(Q) C H — > H be a self-adjoint 
operator on H. For a G {—,0,+} we denote by S a (Q) the set of all nontrivial finite dimensional 
subspaces V C dom(Q) such that Q\ v is negative definite, identically zero, and positive definite, 
respectively. Define! 

fn \ / sup{dimF : V G S a (Q)}, if<S a (Q)^0 , ^ x 

(2.6) a a (Q) = | o, ifS a (Q) = + aG {-'°'+l- 

We introduce the signature of the self-adjoint operator Q by 

(2.7) a(Q):=(a-(Q),a (Q),a+(Q)). 
Recall now the following fundamental observation. 

Proposition 2.5 (Sylvester's Law of Inertia). The signature a is invariant under coordinate 
changes, i.e., o~(Q) = a(S*QS) for any bounded invertible coordinate change S. 

Proof. Invariance follows immediately since the change of variables S is invertible and 

V G S a (Q) implies that S 1-1 ^) G S a (S*QS) and dim^ = dimS _1 (F); 

U€S a (S*QS) implies that S(U) G S a (Q) and dimU = dim 5(17). 



□ 



In the sequel, we will need the following proposition: 



Proposition 2.6. Assume Q : dom(Q) C H — > H is a self-adjoint operator. The the following 
assertions hold true: 

(i) If B is a self-adjoint, invertible, s x s matrix and R : dom(Q) ®C S — > Hc?>C s is defined by 

R= w) b)' then 

(2.8) a a {R) = a a (Q)+a a (B), aG{-,0,+}. 

(ii) If a a (Q), a G {—,0,+}, is finite then it is equal to the dimensions of the stable, zero, and 
unstable eigenspaces of Q, respectively; 

Proof, (i) First we note that our assertion is true if B is a nonzero real number. By induction, 
we conclude that (|2.8p is also true provided the matrix B is a diagonal matrix with nonzero real 
entries. 



Note that in the case of oo, we require Q, not the quadratic form (•, Q-), to be identically zero on V so that ao{Q) 
is the dimension of the kernel of Q. 
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Next, we look at the general case when B is a general self-adjoint matrix. Then, there exists 
Pb an orthogonal matrix and D a diagonal matrix such that B = P^DPb- Since any orthogonal 

matrix is unitary we have that Pg = P B , (^q^ P 1S mver ^i rj le and 



R 



j h or (A 0\ I U 

P B D P B 



Assertion (i) follows from Proposition l2.5l and the observation that (12, 8p is true for diagonal matrices 
B with nonzero real entries. 

We note that assertion (ii) is trivial for a = 0. Assume that a+(Q) is finite. Then, there exists 
a subspace H, a self-adjoint, negative-definite operator Q : dom(Q) C H — > H and D a diagonal 
s x s matrix with positive real entries, where s is equal to the dimension of the unstable space of 

Q, so that Q = Since <r+(Q) = and cx+(D) = s, our conclusion follows immediately 



x Z>, 

from (i). □ 
Assume now 



(H2) Near (X*,u*), there is a C* 1 family (X, c)(u) of critical points of (f!T2|) . 

This follows, for example, by the Implicit Function Theorem, applied to (|2.4p if ker-A^* is trivial 
and also A w * is uniformly invertible. Another situation in which (|H2|) follows, of particular interest 
in the applications considered here, is when ker A w * is due to group symmetries of the problem 
leaving both energy £ and constraints C invariant, and also A L0 *\ is uniformly invertible, 

I (ker 

i.e., A w * has a spectral gap. In general, it may be obtained as part of an existence theory for 
Euler-Lagrange equation (|2.4p carried out by other means. We note that (|H2|) in any case implies 
also 

(2.9) ker A.* 1 VC(X*), 

(by Fredholm alternative), as holds automatically when ker A Wt is generated by group invariances. 
Assume also 

(H3) The matrix (dc) 'du))(u)*) is invertible. 

Differentiating the Euler-Lagrange relation (|2.4p with respect to u, we obtain the key relation 

(2.10) Au*{dX/dojj){u*) = -VCjiX*) for all j = 1, . . . ,r 
and its consequence 

(2.11) (dc k /du 3 )(u*) = -((dX/dLU 3 )(LU*),A w *(dX/dLu k )(LU*)) for all j, k = 1, . . . , r, 

from which we see that (dc/duj)(uj*) is a symmetric form. Our main result is the following identity 
relating the signatures of -A w *|vC(x*)- L ) A^*, and (dc/du)(uj*). 

Theorem 2.7. Assuming ([Hi]) - (jH3j) . 

(2.12) ^(A;*) = ^(^IvC(JC*)-l) +ff(-9c/a.;)(w*)). 

Proof. To prove the theorem we use Sylvester's Law. First, we note that from (|H3p . (|2.10p and 
(|2TTT]> it follows that 

(2.13) H = VC(X*) ± Sp&n{(dX/dLUj)(Lu*) : j = 1, . . . ,r}. 
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Next, we denote H-yCtx*) 1 - an< ^ ^-Span{(dX/dwj)(uj*y.j=i,...,r} the orthogonal projections onto VC(X*) ± 
and Sp&n{(dX/dojj)(uj*) : j = 1, . . . , r}, respectively. Using ([2J3]) we infer that S : VC(X*) ± x 
Span{(dX/duj)(ui*) : j = 1, . . . , r} — >■ H defined by 

(2.14) 5 =h 1 + h 2 



is bounded, invertible with bounded inverse. Moreover, from (|2.10p and (|2.1ip we obtain that 

(2-15) n vc -( X *)X^4a;*nsp a n{(ax/aa)-,)(a)*):j=l,...,r-} = °- 

Since A^* and the orthogonal projections are self-adjoint operators, passing to the adjoint in (|2.15p 
yields 

(2-16) n Span {( a x/awj)(a)*):j=l,...,r}^*n vC ( X «)x =0. 

Finally, from (|2.14p - (|2.16p we conclude that 
(2.17) 

g*A m g = ( ~Rvc(x*)- s -Aw*\vc{x*)- s - 



^ Span {( a x/^w ^ )(a)*):J•=l ) ... ) ^}A)*|spa^l{(aX/^u; j )(w*):j=l,■••.»•} 
Next, we note that from (|2.10p and (|2.1ip it follows that 

r 

(2.18) n Span{(ax/a ^. )(w * ):j=lv .. jr} A^/i = - Yl (dck/dujj)(u;*)(h, (dX/du k )(u*))(dX/duj)(u*), 

j,k=i 

which implies that 

(2-19) ^rns pan {(ax/8w 3 0(w*):j=l,...,r}A>;*|Span{(SX/aw i )(u;*):j=l,... ) r}) = CF^-{dc/ dio){bJ*)^ . 

Since a(Il^r C ^ x ,-j± A^* | vc(x*) x ) = CJ (^w* I vc(x*) x )> the theorem follows shortly from Proposi- 
tion E^i), (ESP and ([2TT9D . □ 



2.3. Variational stability. Denoting ^ := dimker A u *, we strengthen ()H2|) to 

(H2') Near there is a C 1 family (X,c)(u,s) of critical points of $T2\), s G R*. 

With (|H3p . this implies that, for fixed c, the set of minimizers near X* forms a C 1 i'-dimensional 
family 

(2.20) {X(w(c),s) : s G R € }. 

Differentiating (|2.20p with respect to s, noting that C(X(w(c), s)) = constant, and recalling (|2.4p . 
we find that 

(2.21) £(X(a;(c), s)) = constant along the family of nearby critical points. 

Definition 2.8. A critical point X* of ([2.1 p is locally strictly orbitally stable in H C H if, for 
l£ H in a sufficiently small H-neighborhood of X* , S(X) > £(X*), with equality only if X lies 
on the set of nearby critical points, and orbitally stable in H C H if, for X £ H in a sufficiently 
small H-neighborhood of X* , £(X) > £{X*). 

We add the further hypothesis, denoting H-distance from v £ H to a set S C H by d jj(t>, S*): 
(H4) 

) = for # > 0, all v £ VC ± (X*), djj^ker^i^^) < K(v, A^v) H . 



Corollary 2.9. Assuming ([HI]) . f|H2'j) . (fH3|) . and ([H4]) . a croftcai pomi X* o/ (pTj) is strictly 
orbitally stable in Diff(£) if and only if 

(2.22) a-(A u *) = a-(^(-dc/duj)(uj*)J 

and dimker A Uit \\/c- L (x *) = ^> i- e -> family (jH2|) is transversal, and orbitally stable in Diff(£) only 
if 

(2.23) MA,*) = a-({-dc/duj){u*)), 
where uj* is the (unique) Lagrange multiplier associated with X* . 

Proof. Immediate, from (j2. 12|) . (j2.21[) . and (|H4p . together with a second-order Taylor expansion 
like that of Proposition 12.21 along perpendiculars to the manifold (jH2]) where C = c*. □ 

Remark 2.10. In the PDE context, Assumption (jH4|) effectively limits the regularity of the space 
H. More generally, the H-norm cannot be stronger than any norm under which A Wt is bounded, 
or (|H4|) would fail. If A u , is bounded on H (in particular, if Diff (£) = H), then H is determined. 
In applications, ()H4p typically follows by coercivity of A Ut , in the form of a Garding inequality, 
together with the presence of a spectral gap. For all of the (PDE) examples of this paper, there is 
a unique space H on which A Wt is both bounded and coercive in the strong sense of (|H4p . 



3. Time-evolutionary stability 
Corollary 12.91 yields stability results for time-evolutionary systems 
(3.1) (d/dt)X = F(X) 

with variational structure, in particular, for energy conserving (e.g., generalized Hamiltonian) , and 
energy dissipating (e.g., quasi-gradient type) systems, as we now describe. 

For this section, it will be convenient to introduce a graded sequence of Hilbert spaces Hi C 
H C H2, with associated inner products (■,-}h^ ) (v)h of varying strength, and gradients V^, 
satsifying: 

(Bl) For X G Hi, Vh 2 and V^ 2 of £ and C are well-defined. 

In particular, this implies that Hi C Diff(£). 

3.1. Generalized Hamiltonian systems with constraints. We first recall briefly the context 
considered in |GSS] of a generalized Hamiltonian system 

(d/dt)X = J(X)V U2 £(X), t>0, 

[6 ' x(o) e H, 

J = — J* , which has additional conserved quantities 

(3.3) C(X(t)) = c G W, i.e., (J(X(t))VC 3 (X(t)),V£(X(t))) = 0, j = 1, . . . ,r, t>0 
besides the Hamiltonian £(X(t)), which is always conserved, by skew-symmetric structure 

(3.4) (V£(X(t)),J(X(t))V£(X(t))) =0. 

Here J = —J*, is assumed to be a closed, densely defined linear operator with X-independent 
domain. 

For such systems, constraints may be divided into two rather different types. The first type, 
familiar from the classical mechanical case, are constraints Cj : H — > R for which 

(3.5) gj {X) := J(X)VCj(X) ^0, for all j = l,...,r. 



These are associated with group symmetries Gj(0, •) leaving the Hamiltonian 8, and thus the flow 
of ([3.2p . invariant, generated by the gj, through 

(3.6) (d/dB)G j (0;X)=g j (G j (9;X)), G j (0;X)=X, j = l,...,r. 

When J is one-to-one0 these are the only type of conserved quantities that can occur, and ev- 
ery group invariance determines a conserved quantity, following Noether's Theorem, through the 
relation VCj(X) := J^g^X), where gj (X) := G'j(0;X), j = l,...,r. 

When J is not one-to-one, however, as for example in the KdV and Boussinesq cases considered in 
|GSSllZ2|IBJK| . for which J = d x , conserved quantities may or may not correspond to symmetries. 
For related discussion, see |Z2j . In this case, we identify a second class of constraints, closer to the 
themes of this paper, for which 

(3.7) JVCj = for all j = 1, . . . , r. 

These are not associated to a nontrivial group symmetry, since gj = JVCj, j = 1, . . . ,r generates 
the identity. 

For example, in the case considered in jGSSj of KdV on the whole line, r = 1, and the constraint 
is associated to a group invariance G(s) consisting of spatial translation by distance s. For the 
generalized KdV case considered in |BJK| on bounded periodic domains, r = 2, only one of the two 
constraints corresponding to translation invariance and the other corresponding to conservation of 
mass (integral over one period). 

Denote the (possibly trivial) symmetries generated by constraints Cj through (|3.6h as Gj(9;X), 
j = r. 

Following |GSS| . we assume that (|3.2p has a well-defined solution operator in H conserving the 
formal invariants C and £: 

(B2) Equation (|3.2p admits a (weak) solution X(t) £ H continuing so long 

as ||X(t)|| jj < +oo, with C(X(t)) and £(X(t)) identically constant. 

Definition 3.1. A generalized traveling- wave solution of (I3.2D is a solution whose trajectory is 
contained in the orbit of its initial point V(0) under the action of the invariance group generated 
by {Gj : j = l,...,r}. 

Lemma 3.2. A solution emanating from X(0) 6 Diff(£) is a generalized traveling wave if and 
only if X(0) is a critical point of the constrained minimization problem (I2.ip . in which case X(0) 
is a stationary solution of the shifted equation (d/dt)X = J(X)V£ UJ * (X), with associated linearized 
operator 

(3.8) (V(0)) := <LF U . (V(0)) = J(X(0))A u * . 

Proof. By group invariance of (|3.2p . it is equivalent that X'(0) = J(X(0))\/£ UJ * (X(0)) lie in the tan- 
gent space Span{ J(X(0))VCj(X(0)) : j = 1, . . . , r} of the invariance group. Since ker J(X(0)) C 
Span{VCj(V(0)) : j = 1, . . . ,r}, this is equivalent to (|2,4p . which in turn is equivalent to station- 
arity as a solution of the shifted equation. Writing, equivalently, T = JV£ UJ * , differentiating in X, 
and using (|2.4p . we obtain ()3.8|) . □ 

Definition 3.3. A generalized traveling wave X(t) € Hi is Hi — > H orbitally stable if for every 
e > 0, there is 5 > such that, for V(0) within 5 distance in Hi of V(0), Y(t) exists for all t > 
and remains within e distance in H from the orbit of X(0) under the group of invariances generated 
by {Gj:j = l,...,r}. 

Note: by skew-symmetry, J is one-to-one if it is onto, as sometimes assumed |GSS) . We remark that J is also 
Fredholm in the typical situation where we have such non-Noetherian constraints, for example in the case of NLS 
(J = i) or gKdV (J = d x ) on bounded periodic domains, in which case J is one-to-one if and only if it is onto. 
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Given a generalized traveling wave solution X(t), we assume further: 
(CI) kerL„*(X(0)) C Span{ J(X(0))VCj(X(Q)) : j = 1, . . . ,r} (transversality) , 

that is, that all neutral directions are accounted for by group symmetry. 

Proposition 3.4. Assuming (|Hl|) . (|H2' p .pI3 |> - ((H4|) . pT ]l -(|B2 |) , (|CT|) . a generalized traveling-wave 
solution X(t) of (j3.2[) is H — >■ H orbitally stable if X(0) G Hi is a strict local orbital minimizer of 
(|2,ip . or, equivalently, the signature of (dc/duj)(uj*) satisfies condition (|2.22p . and Diff(£) = H. 

Proof. By (|H1|) . the level sets of C smoothly foliate i?(X(0),5), whence it is sufficient to consider 
perturbations such that C(Y(0)) = C(X(0)). As £(Y(t)) — £(X(t)) is constant, and, by strict 
orbital minimality, majorizes the || • ||jj distance of Y(t) from the manifold of equilibrium solutions 
through X(0), which, by (ICip . is the group orbit of X (0)0 we obtain global existence by (|Bip and 
continuation, and e-boundedness by the assumed continuity of £ with respect to || • ||h- D 

Remark 3.5. The notion of orbital stability defined above does not include phase-asymptotic sta- 
bility; that is, it does not assert that Y(t) stays close to X(t). It is not difficult to see by specific 
example (for instance, the different-speed traveling-waves of [GSS]), that this is all that we can say. 

Remark 3.6. Condition (jCip could be weakened to the assumptions that for c near c* the manifold 
of critical points (]2.20p extends globally as the orbit under the group symmetries generated by 
{Gj : j = 1, . . . , r} of a (for example, compact) manifold of strict local orbital minimizers of (|2.ip 
on which A w has a uniform spectral gap, a possibility that does not seem up to now to have been 
discussed. 

Remark 3.7. The "one-way" stability result of Proposition 13.41 recovers the result of [GSS for 
solitary waves of KdV; for J onto (one-to-one), a converse is given in |GSSj for r = 1. In general, 
time-evolutionary stability does not imply constrained variational stability but only certain parity 
information on the number of unstable eigenvalues. For complete stability information, a more 
detailed analysis must be carried out, such as the Krein signature analysis in |B JK] . 

3.2. Quasi-gradient systems with constraints. We now turn to the situation of our main 
interest, of quasi-gradient systems 

(3.9) (d/dt)X = -M(X)V£(X), -(M + M*) > 0, 

with constants of motion C(X(t)) = cq 6 W for which 

(3.10) (d/dt)£(X(t)) = {V£(X(t)),M((X(t))V£(X(t))) < 0, 
guaranteed by 

(3.11) MVCj = for all j = 1, . . . , r, 

similarly as in (|3.7p . Recall, for M semidefinite, ker M = ker M* C ker ^(M+M*), by consideration 
of the semi- definite quadratic form \{M + M*), so that (|3"TTD implies also (VCj,MV£) = for 
all j = 1, . . . , r, giving (d/dt)C(X(t)) = 0. 

For the Keller-Segel (chemotaxis) model considered in |PSH IPS2j . r = 1, with the single con- 
straint corresponding to conservation of mass (total population). For the equations of viscoelas- 
ticity with strain-gradient effects on bounded periodic domains, considered in |OZH IY1 IBYZj . 
r = 2d = 2,4,6, where d = 1,2,3 is the dimension of allowed deformation directions, with the 2d 
constraints corresponding to conservation of mass in deformation velocity and gradient coordinates. 



Here, we are using group invariance to conclude from local majorization a global property as well; this, precisely, 
is the step where we require (|Clf) . 
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Lemma 3.8. X(t) = X{0) £ Diff(£) is a stationary point of (|3.9p if and only if it is a critical point 
of the constrained minimization problem (|2.ip , in which case the associated linearized operator is 

(3.12) L u * (X(0)) := dF w . (X(0)) = -M(X(0))A u * . 

Proof. Essentially identical to the proof of Lemma 13.21 □ 

For the applications we have in mind, we find it necessary to substitute for (|B2[) the milder 
assumption: 

(B2') For data X(0) £ Hi near an equilibrium solution X*, (|3.ip admits a strong 

solution X(t) £ Hi, with (d/dt)X(t) = T{X(t)) £ H 2 . Moreover, this solution 
continues so long as X(t) remains sufficiently close in H to an equilibrium. 

Remark 3.9. Condition (|B2'P is satisfied for all of the example quasigradient systems considered 
in this paperQ whereas (|B2p is satisfied only for the Cahn-Hilliard equations considered in Section 
PI 

3.2.1. Orbital stability. Assume analogously to (|C1|) : 

(CI') kerL^* is generated by group symmetries of (|3.9p . 

Definition 3.10. A stationary solution X(t) = X(0) £ Hi of (|3.9p is orbitally stable if for every 
e > 0, there is 5 > such that, for Y(0) within 5 distance in Hi of X(0), Y(t) exists for ail i > and 
remains within e distance in H of the manifold (|2.20p of stationary solutions (critical points of (|2.ip ) 
near X(0). It is asymptotically orbitally stable if it is orbitally stable and, for ||5^(0) — X(0)||hi 
sufficiently small, Y(t) converges in H to the manifold (|2.20p of stationary solutions (critical points 
of flU}) near X(Q) having the fixed value c = C(Y(0)). 

Proposition 3.11. Assuming ([HI]) . pLT]) . (fH3|) - ([Hi]) . (|EI j) -(|E2* j) . (fCT]) . a stationary solution 
X(t) = X(0) of (|3.9|) is Hi —7- H orbitally stable if X(0) £ Hi is a strict local orbital minimizer 
of (|2.1h . or, equivalently, the signature of (dc/dco)(uj*) satisfies condition (j2.22j) . 

Proof. Identical to that of Proposition 13.41 □ 

3.2.2. Asymptotic orbital stability: finite- dimensional case. Restricting for simplicity to finite di- 
mensions, take without loss of generality H 2 = H = Hi, by equivalence of finite-dimensional 
norms. 

Now, make the additional assumption: 
(Dl) Constant-energy solutions 8{X{t)) = S(X(Q)) of are stationary: X(t) = X(0). 
Lemma 3.12. A sufficient condition for ()Dip is M = M* , or, more generally,. 
(Dl') kerM = ker ±(M + M*), 

Proof. Using (|3. lOj) we infer 

(d/dt)£(X(t))\t=o = -(V£(X(0)),M(X(p))V£(X(p))) 

= -l(VS(X(0)), (M + M*)(X(0))V£(X(0))) 

vanishes if and only if V£(X(0)) lies in ker(±(M + M*)(X(0)) \ = kerM(X(0)), whence 

T(X(0)) = -M(X(0))V£(X(0)) = 0. 
Noting, finally, that JF^* = J 7 , we are done. □ 



^ The continuation property follows by nonlinear damping (energy) estimates as introducted in [Z4I IZ1| . which 
show that higher Sobolev norms are exponentially slaved to lower ones.. 
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Remark 3.13. As (jDl'P shows, (|D1|) is typical of gradient systems. That is, gradient systems 
typically do not possess nontrivial constant-energy solutions such as time-periodic solutions or 
generalized traveling- waves, in contrast with the situation of the Hamiltonian case. 

Proposition 3.14. In finite dimensions, assuming (|Hip . (|H2'|) . (|H3p . (|Cl'p . and (|Dip . a station- 
ary solution X(t) = X(0) of (|3.9p is H — > H asymptotically orbitally stable if X(0) G H is a strict 
local orbital minimizer of (|2.ip . or, equivalently, the signature of (dc/duj)(co*) satisfies condition 
([232]) . 

Proof. First note that (|H4p and (|Bip - (|B2'P hold always in finite dimensions. By orbital stability, 
Y(t) remains in a compact neighborhood of the orbit of A(0) under the group symmetries of (|3,9D . 
whence, transporting back by group symmetry to a neighborhood of A(0), we obtain an w-limit 
set which is invariant up to group symmetry under the flow of (I3.9p . Moreover, by nonincreasing 
of £(Y(t)), and local lower-boundedness of £, we find that £ must be constant on the w-limit set. 
Combining these properties, we find that points of the cu-limit set must lie on constant-energy 
solutions, which, by (jDlh . are stationary points of (13. 9p . Noting that, up to group symmetry, Y(t) 
approaches its w-limit set by definition, we are done. □ 

Definition 3.15. A stationary solution X(t) = X(0) G Hi of (|3.9p is Hi — > H stable if for every 
e > 0, there is 5 > such that, for Y(0) within 5 distance in Hi of -A(O), Y(t) exists for all t > 
and remains within e distance in H of A(0). It is Hi — > H phase-asymptotically orbitally stable if 
it is stable and, for ||y(0) — A(0)||hi ^ <^> Y{t) converges time-asymptotically in || • ||h to a point 
X^(Y(0)) on the manifold ()2.20p of stationary solutions (critical points of (|2.ip ) near X(0). 

Proposition 3.16. In finite dimensions, assuming (|HT|) . plT]) . (jH3]) . (|DT]) . a stationary solution 
X(t) = X(0) 6 H of (j3.9[) is H — > H stable only if X(0) is a nonstrict local minimizer of (|2.1D . 
or, equivalently, the signature of (dc/duj)(u*) satisfies condition (|2.23|) . It is phase- asymptotically 
orbitally stable, indeed, time- exponentially convergent, if X(0) is a strict local minimizer of (I2.ip . 
or, equivalently, the signature of (dc/duj)(uj*) satisfies condition (12.221) . 

Proof. Again, note that (IH4h and (IBip - (IB2'P hold always in finite dimensions, (i) If A(0) is not a 
minimizer, then there is a nearby point Y(0) withC(Y(0)) =C(X(0)) = c* and £(Y(0)) <£(X(0)). 
As C(Y(t)) = C(y(0)), the only stationary points to which Y(t) could converge are on the manifold 
(I2.20p with \s\ small and c = c*, along which £ = £(X(0)). But, then, by continuity of £ in H, 
£(Y(t)) would converge to £ (A(0)), a contradiction. This establishes the first claim. 

(ii) By (IHip . the level sets of C smoothly foliate B(X(0),5), whence it is sufficient to consider 
perturbations such that C(Y(0)) = C(X(0)). If X(0) is astrict local minimizer, then, in B(X(0), e), 
£(Y(t)) majorizes the || • ||-distance from Y(t) to the manifold of stationary solutions through X(0). 
Moreover, 

(d/dt)£(Y(t)) = (V£(Y(t)),M(Y(t))V£(Y(t))} < -||n vc(x(0)) x V£(Y(t))\\ 2 < -(SQr(t)-S(X% 

where ^-^cixio)) 1 - denotes orthogonal projection onto VC(X(0)) ± , and X^ denotes the orthogonal 
projection of Y(0) onto the set C(X) = cq. Here, the final inequality follows by the assumed 
variational stability. Thus, £(Y(t)) — £(X(0)) decays time-exponentially so long as Y(t) remains 
in ^(A^O^e), as does therefore ||n vc .( X ( ^± V£(Y(t))\\. Thus, 

j t ((dX/ds)\ s=0 , c *,Y(t) -y(0)) < K\\M(Y(t))V£(Y(t))\\ < Ke^\ for some K,rj > 0, 

and so ((dX/ds)\ s= o tC * , Y(t) — Y(0)) converges time-exponentially, which, together with the al- 
ready established convergence to the manifold of equilibrium solutions, gives stability and phase- 
asymptotic orbital stability, with convergence at exponential rate. □ 
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3.2.3. Asymptotic orbital stability: infinite- dimensional case. Under suitable assumptions guar- 
anteeing compactness of the flow of (|3.9p . the Lyapunov argument of Proposition 13.141 may be 
generalized to the infinite-dimensional case, yielding orbital stability with no rate; see |Wr| . This 
applies in particular to all of the examples given below, which are (or are equivalent to) systems 
of parabolic type on bounded domains. However, note that this requires a spectral gap of A^*, 
under which we would expect (assuming a correspondingly nice M, as also holds for our examples), 
rather, exponential stability, whereas the exponential stability result of Proposition 13. 161 first, does 
not generalize in simple fashion to infinite dimensions, due to the presence of multiple norms, and, 
second, requires the restrictive assumption (jDl'l) . which does not apply the important example of 
viscoelasticity, below. 

Moreover, there are interesting situations in which the essential spectrum of A^* and cLF, extend 
to the imaginary axis; for example, stability of solitary waves in viscoelastic, Cahn-Hilliard, or 
chemotaxis systems, as studied respectively in |Z3 1IHoZ] . |Hol|,lHo3] . and |PS2] . In such situations, 
it has proven useful to separate the questions of spectral and nonlinear stability; see [Zl] for a 
general discussion of this strategy. In particular, for each of the above-mentioned settings, it can 
be shown that to show nonlinear orbital stability, with time-algebraic rates of decay, it is enough to 
establish transversality plus strict stability of point spectrum, whereupon essential spectrum and 
nonlinear stability may be shown by separate analysis @ 

For both of these reasons, we restrict ourselves in the infinite-dimensional setting to the simpler 
treatment of stability of point spectrum of dJ~, from which the relevant nonlinear stability result 
may in most settings (in particular, all of those mentioned in this paper) may be readily deduced. 

Definition 3.17. A stationary solution X(t) = X(0) G Hi of (j3.9|) has orbitally stable point 
spectrum if (i) the dimension of the zero eigenspace of dJ-(X(0)) = M(X(0))A UJ * is equal to the 
dimension dim ker A* + dime of the manifold (|2.20p of stationary solutions (transversality), and 
(ii) the nonzero point spectrum of dJ-(X(0)) has strictly negative real part. 

Definition 3.18. A stationary solution X(t) = X(0) G Hi of (|3.9p is Hi — > H-linearly asymptoti- 
cally orbitally stable if, for Y(0) G Hi, the solution Y(t) of (d/dt)Y = M(X(0))A w *Y converges in 
H to the tangent manifold of the manifold ()2.20p of stationary solutions (critical points of (|2.ip ). 

Define now the following local linear analog of the nonlinear condition (IPip : 

(El) The center subspace of M(X(0)) A;* consists of ker M(X(Q))Au* ; 

that is, is the only pure imaginary eigenvalue of M(X(0))A u * , and associated eigenvectors are 
genuine. 

Lemma 3.19. Sufficient conditions for (|Eip are (jPl'P (e.g., M = M* ) and A u * | vc*(X(o)) i - — 0- 
(jDl'P by itself implies that there are no nonzero imaginary eigenvalues of M(X(0))A Uti . 

Proof. Let M(X(0))A^Y = ifiY, /*£l Then, dropping uj* and X(0), (AY, MAY) = ifi(Y, AY), 
hence 

(AY, (M + M*)AY) = 2Re(AY, MAY) = 0, 
and, by semidefiniteness of (M + M*), we conclude that (M + M*)AY = 0. By (jDl'p . this implies 
MAY = 0. This shows that is the only pure imaginary eigenvalue of MA. 

Now, suppose that (MA) 2 Y = 0. Then, by §DV\ . M*AMAY = 0, hence (MAY,A(MAY)) = 0. 
As MAY G VC(A(0)) ± , by ([3TTT]) , and ^| V c(x(o))-l > by assumption, this implies that MAY G 
ker A, hence (AY, (M + M*)AY) = and so AY G ker(M + M*) = ker M and MAY = as before. 
This shows that there are no generalized eigenvectors in the zero eigenspace of MA, completing 
the proof. □ 



In the periodic parabolic examples studied below, nonlinear stability follows easily from spectral stability by 
general analytic semigroup theory [He I [S] . 
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Remark 3.20. Transversality condition (i) of Definition 13.171 concerns well-behavedness of the ex- 
istence problem. Under (|E1|) . it is guaranteed by our hypotheses (|Hip - (|H4|) . through (|2.10p . and 
the associated relation ker M(X(0))A U)t = hsr A^* © Sp&n{dX/ du)j)(u)*) : j = 1, . . . ,r}. 

Along with (|Eip . we identify the following mild properties of the linearized solution theory H 
(E2) If X(0) is asymptotically orbitally stable, then solutions of (d/dt)Y = M(X(0))A UJ *Y 
with Y(0) in the constrained unstable subspace of A u „ exist for all time and converge in H to the 
center subspace of M(X(0))A UJ * . 

(E3) Au„ is bounded on H as a quadratic form. 

Theorem 3.21. Assuming (jHl|) ; ([H^ . (fH3|) - (fH4]) . ([BTjl - lfTW]) . and (jEl|) ; a stationary solution 
X{t) = X(0) G H of (|3.9p has orbitally stable point spectrum if and is H — > H linearly orbitally 
asymptotically stable only if X(0) is a strict local minimizer of (|2. 1 j) . or, equivalently, the signature 
of (dc/duj)(uj*) satisfies condition (|2.22p . Assuming also (|E2p - (|E3p . X(0) has strictly orbitally 
stable point spectrum if and only if X(0) is strictly orbitally stable, or (|2.22p . 

Proof (<*=) (i) Dropping u* and X(0), MAY = is equivalent to AY = E^=i «i v Ci(^(°))> for 
some a G C r . Since A^-(u*) = -VCj(X(0)) for all j = 1, . . . , r, we thus have 

3=1 J 

or Y G Spaa{SJT/ajj-(w*) : j = 1, . . . ,r} © Span{VC j (X(0)) : j = 1, . . . ,r} as claimed. 

(ii) By (jElj) . it is sufficient to show that ReA < whenever -M(X(Q))A UJ *Y = XY. Dropping uj* 
and Z(0)) again, we find from -Re(AY, = ReA (Y, AY) that either AY = 0, hence A = 0, 

or ReA = zffiffiS < 0. 

(=>) On the other hand, suppose X(0) is not strictly orbitally stable. Then, either ^4 W * I vc(x(o)) ± 
has a zero eigenvector not in the tangent subspace (dX/ds) of the family of equilibria described 
in ([230]) . violating transversality, or else £(Y(0)) := (Y(0), A^Yfi)) < for some F(0) G 
VC(A"(0))- L , in which case we see from (d/dt)(Y(t), A u *Y(t)) < that £(F(t)) < for all t > 
for the solution Y of the linearized equations with initial data Y(0). Thus, Y(t) cannot approach 
the zero-£ set tangent to the manifold of nearby minimizers as t — > oo, violating linear asymptotic 
orbital stability. Assuming also (|E2p . so that we have transversality by Remark I3.20| we find that 
there are strictly unstable spectra of M(X(0))A ult , or else (|E1|> — (|E3[) would yield convergence to 
the manifold of nearby minimizers, a contradiction. □ 

Remark 3.22. It would be interesting to establish a relation between the number of constrained un- 
stable eigenvalues of A u * and the number of unstable eigenvalues n^(M(X(0))A u) * ) of the operator 
—M(X(0))A ul * , similar to those obtained for A u * and J(X(0))A w * in the Hamiltonian case (see, e.g., 
|GSSl IBJKj ) . These are equal (under appropriate assumptions guaranteeing that eigenvectors re- 
main in appropriate spaces under the action of M(X(0)) ±1 / 2 ) when M(X(0)) is symmetric positive 
definite, by the similarity transformation M(X(0))A w * -> M(X(0))~ 1 / 2 M(A(0)) J 4 aJ »M(X(0)) 1 / 2 = 
M(X(0)) 1 / 2 A w .M(X(0)) 1 / 2 . More generally, they are equal if M(X(Q)) + M(X(0))* > 0, by homo- 
topy to the symmetric case, together with the observation that, by (jDl'P and Lemma 13.191 zero is 
the only possible imaginary eigenvalue of M(X(0))A UJ * , and, since M(X(0))A LJ *M(X(0))A UJ *Y = 
implies A^M^X^A^Y = and thus fa(^*Y,M(X(Q))A u .Y) = 0, hence Y G ker A,*, from 



This is satisfied for all of the examples of this paper, not only on the unstable subspace, but all of H as expected. 
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which observations we may deduce that no eigenvalues cross the imaginary axis throughout the 
homotopy. 

We conjecture that, more generally, under assumptions (|Dr|) and (|Hip . (jH2'|h(|H3|) . 

(3.13) n-(M(X(0))A u *) = <x-(A,*lvc^(x(o))) = MA.*) - *-(dc/dw)(u*)), 

where (jH3j) is needed only for the second equality. In the finite-dimensional, diagonalizable, 
case, this follows easily from the observation that M {X{Q))A W *Y = fiY, Re// ^ implies that 
Y e VC(X(0)) x , Y kerA,*, and {A W *Y, (M(X(0)) + M(X(0))*)A LU *Y) > 0, and, moreover, 
Re(A UJ *Y, M(X(0))Au*Y) = Re/x(Y, A^Y), so that {Y,Au*Y) ^ 0, whence n^(M(X(0))A lJ *) < 

(X*|vC(X(o))^) ■ and n + (M(X(0))A^) < <7+(^*| vc?(x(0)) ±J . Recalling that dHU) , dHU) (via 

(gSD ) imply n (Mpf(0))A/0 = dimker M(X(0))A,* = 0o(A,*) + dimker M(X(0)), we thus ob- 
tain the result. Galerkin approximation (with well-chosen subspaces including ker(M(X(0))A aj * )) 
then yields the result in the case that M(X(0))A LJ * has a spectral gap and finitely many stable 
eigenvalues. 

4. Examples: applications to co-periodic stability of periodic waves 
4.1. Cahn Hilliard systems. First, we consider the Cahn-Hilliard system 
(4.1) u t = dl(VW{u)-d 2 x u) 

where u G R m and W : W 71 — >■ M is a smooth function. This model can be expressed in the form 
(I3T9D with 



(4.2) £(u) := J (\\ u x\ 2 + W{ufjdx and M := -fig > 

(4.3) V L2 £:(n) = VVF(u) - ^u. 

Note that scaling x, t, and W, we can restrict to the unit interval without loss of generalization. We 
note that the energy is continuously differentiable on H^ cr [0, 1] and L 2 -gradients are continuous on 
Hp Cr [0, 1]. The system fits our framework so that we can derive stability and instability results that 
recover and substantially extend the spectral stability conclusions of |Hol] for spatially periodic 
solutions of (|4.ip . We define the constraint function by 

C(u) = j u{x)dx. 
Jo 

To find solutions of (|2.4p for the Cahn-Hilliard case, we need to solve the nonlinear oscillator 
problem 

(4.4) -n« + VlfW = -w*, 

with uj* E M. m . Whenever the potential W is not convex, ()4.4p possesses families of periodic 
solutions which are in fact critical points of (]4.ip . For the linearization at such periodic solutions, 
we find simply the second derivative of the energy since the constraint function C is linear: 

(4.5) A.* = V| 2 £(u*) - lo*V 2 l2 C(u*) = -d 2 x + V 2 W(u*). 

In the scalar case, m = 1, consider periodic patterns with minimal period 1, first. These possess 
precisely 2 sign changes, so that by Sturm-Liouville theory, the translational eigenfunction u* is 
either the second or third eigenvalue of A u *. More precisely, one can decompose the spatially 
periodic linearized problem on [0, 1] into a problem on even and odd functions, which, equivalently, 
satisfy the linearized problem on [0,1/2] with Neumann and Dirichlet conditions, respectively. 
Since u* has a sign on [0,1/2], the linearization is semi-definite on the Dirichlet subspace. The 
Neumann problem can have cr_ = 1 or = 2. It is well-known that the index in the Neumann 
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case depends on the period-amplitude relation. In fact, periodic orbits patterns in the scalar case 
come in families parameterized by the maximum u max of the solution, once we allow the spatial 
period L to vary. Whenever the period is increasing with amplitude, we find a = — 1, decreasing 
period corresponds to a = —2. We refer the reader to |Korl ISchaj and references therein for proofs 
in the analogous Dirichlet case. Alternatively, this could be (but to our knowledge has not been) 
deduced by a co-periodic stability index computation like those in |OZll IHolj in the conservation 
law case, which would yield that the Morse index be even or odd accordingly as du max /dL (related 
by a nonvanishing factor to the derivative of the Evans function at A = 0) is less than or greater 
than zero. In particular, we find that periodic patterns are co-periodically stable if and only if 
du mgbX /dL > and dc/du) > 0. 

In the system case, we obtain a partial analog, namely, a necessary and sufficient condition for 
co-periodic stability of periodic solutions of (|4.ip is that the number of positive eigenvalues of 
(dc/dco)(co*) is equal to the number of negative eigenvalues of A^* = —d 2 + V 2 W(u*), considered 
as an operator on all of L 2 [0, l] per (i.e., with no zero-mass restriction). This result appears to be 
new, despite the considerable attention to Cahn-Hilliard equations in the pattern-formation and 
phase-transition literature. 

Moreover, for the important case of multiply-periodic solutions in W 71 , which appear (see |Ho2] ) 
not to have been treated so far either in the scalar or system case, we obtain corresponding stability 
conditions by exactly the same argument. Specifically, we find that a necessary and sufficient 
condition for co-periodic stability is that the number of positive eigenvalues of (dc/du;)(uj*) is equal 
to the number of negative eigenvalues of eigenvalues as does A u * = — d 2 + V 2 W(u*), considered as 
an operator on L 2 [0, l] per . These results too appear to be new. 

Finally, we collect the results of this subsection in the following proposition. 

Proposition 4.1. Assume that u* is a periodic or multi periodic pattern of (|4.ip . The fol- 
lowing assertions hold true: 

(i) If m = 1 and x 6 [0, l]p er , the pattern u* is stable if and only if (dc/duj)(u*) > and 
dT I du max 

(ii) If m > 1 or for multi-periodic patterns, u* is stable if and only if —{dc/duj){uj*) has the 
same number of negative eigenvalues as does A^* = —d 2 + V 2 W(u*), considered as an 
operator on L 2 [0, l] per . 

Remark 4.2. Integrating the traveling-wave ODE (|4.4p over one period, we find further that 



where superscript a denotes average over one period of the periodic solution u* , here indexed by 
its mean c = Jq u*(x)dx. This gives the more explicit stability condition 



or a_(d(VW{u*) a )/dc) = a-(A u! *) (case m > 1). 

Remark 4.3. By (|4.6p . (dc/dco)(uj*) simplifies in the long- wavelength limit toward a homoclinic 
solution with endstate u°° as x — > ±oo simply to —W"(u°°), which for m = 1 (by the existence 
criterion that u°° be a saddle point) is scalar and negative, implying co-periodic instability. Stable, 
large- wavelength patterns can be found in the heteroclinic limit, when the periodic solutions limit 
on a heteroclinic loop. In the small-amplitude limit, that is, for almost constant periodic solutions, 
both coperiodic stability and instability can occur. In the simple scalar cubic case, W'(u) = 
L 2 (—u + n 3 ), small-amplitude patterns exist for |c| < 1/V3 for an appropriate length parameter 
L. These small amplitude patterns are stable when |c| < l/\/5 and unstable when |c| > l/\/5; see 



(4.6) 



to* = -VW{u*) 



(4.7) 



d(W'(u*) a )/dc < (case m = 1) 



[GNj. 
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4.1.1. Sturm- Liouville for systems. It is a natural question to ask whether one could obtain in- 
formation about a-{A ul *) also in the system case m > 1 using generalizations of Sturm-Liouville 
theory such as the Maslov index. As described in O |Ej |Uj, there exist such generalizations not 
only for operators of the form A w * , but for essentially arbitrary ordinary differential operators A^* 
of self-adjoint type. However, so far as we know, all of these generalizations involve information 
afforded by evolving planes spanned by d different solutions, rather than the single zero eigen- 
function u* from which we wish to draw conclusions. Indeed, in the absence of additional special 
structure (as described for example in Cor. 15. 2p . we do not know how to deduce such information 
analytically. 

On the other hand, we note that this point of view does suggest an intriguing numerical approach 
to this problem. Namely, borrowing the point of view espoused in (EllUj that Sturm theory amounts 
to the twin properties of monotonicity in Morse index of the operator A^* — Aid on domain [0, x] 
with respect to frequency A and spatial variable x, we find that we may compute the number 
of eigenvalues A < Ao of A^* on domain [0, xq] by counting instead the number of conjugate 
points x G (0,xo) of A^* — Aold, i.e., values for which A^* — Aold has a kernel on [0,x]. This 
in turn may be computed using the same periodic Evans function of Gardner [G] , as described 
for example in [B JNRZ] . that would be used to study behavior in A. Some apparent advantages 
of this approach are that (i) it avoids the computationally difficult large-|A| regime, and (ii) the 
computational expense is no greater in computing values of the Evans function along all of [0, xq] 
than in evaluating at the single value x = xq, which in any case requires integrating the eigenvalue 
ODE over the whole interval. This seems an interesting direction for further investigation. 

4.2. Coupled conservative reaction diffusion equations. We consider the general model 

/ 4 ^ u t = d x (a(u)u x - b(u)v x ^j , 

v t = v xx + Su + g(v), 

where a, b, g are of class C 3 and a(u) > ao > 0, b(u) > bo > 0, u, v 6 R, 5 > 0. We note that this 
system can be expressed in the form ()3.9[) with 

(4.9) S(u,v):=J {F(u)-uv-5- l G(v) + (25)- l \d x v\ 2 ^dx, 
where F"{u) := |M and G'(v) = g(v), and 

,A-, n \ ^ (d x (Uu)d x ) OA ^ n _ . / F'(u)-v 

(4.10) M(u, v):=-^^ Q >*> d j>0, V L2 £(u, v) = \ _ § _ lg{v) 

This general form of a conservation law coupled to a scalar reaction-diffusion equation includes a 
variety of interesting models such as the Keller-Segel equations of chemotaxis and many of their 
variations, reaction-diffusion systems with stoichiometric conservation laws, super-saturation mod- 
els for recurrent precipitation and Liesegang patterns, as well as phase-field models for phase 
separation and super-cooled liquids; see for instance |KSt IKR[ IPS21 |PS3] . 

Again, these systems fit into our general framework of quasi-gradient systems and we can relate 
spectral stability to signatures of symmetric forms. Since the system f|4.8j) is parabolic, spectral 
stability will readily imply nonlinear asymptotic stability up to spatial translations. We thereby 
recover and extend stability computations from [PSlj . 

In this case the constraint function is given by 

C(u, v) = I u(x)dx. 



o 
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Next, we look for solutions of system (|2.4p . We note that in the case at hand this system is 

F'(u)-v = -u*, 

(4.11) 



+ g(vfj -u = 0. 




Since F" > 0, we can solve the first equation of (|4.1ip for u. We then substitute the solution 
u = (p(v — lj*) in the second equation of (I4.1ip . and conclude that the u-component of any solution 
of (I4.1ip satisfies the nonlinear Schrodinger equation 

(4.12) v xx + 8<p(v - oj*) + g(v) = 0. 

The above w*-family of equations is a parametrized nonlinear oscillator, well-known to support, in 
appropriate circumstances, a variety of homoclinic, heteroclinic, and periodic orbits. 
As the constraint function C is linear, we find that 

( a ( u *) 
b(u*) 

In analogy to the case of the Cahn-Hilliard equation, we can give necessary and sufficient condition 
for co-periodic stability of spatially periodic patterns. Denote by w max the maximum of v and write 
L(u max ) for the period of the family of periodic solutions to (|4.12p . 

Proposition 4.4. The necessary and sufficient condition for co-periodic stability of spatially pe- 
riodic patterns of coupled conservative-reaction diffusion equations (|4.8p is {dc/ duj){uj*) > and 
dL/dv max > 0. 

Proof. To start, we consider (u*,v*) a co-periodic solution of (|4.1ip of period 1. We note that all 
we need to show is that a-{A UJ *) = 1. Since (p = (F 1 )" 1 and u* = (f(v* — uj*), we obtain that 

, , 't * *\ _ 1 1 b ( u *) 

[ ' ~ U ' ~ F"(ip(v* - oj*)) ~ F"(u*) ~ 'oJu T )' 

Since v* is the solution of the nonlinear oscillator (|4.12p wg Ccin again conclude that ^4^* '. — 



-$(^ + sV)). 



b(u* ) 

(4.15) A^* 7 — r has only one simple, negative eigenvalue, 

a[u*) 

provided that dL/dv max > 0. 

One can readily check that A^* has the following decomposition: 



(4.16) A UJ *=S*T W *S, where S = M>(«*) u and T L 





-4, 



b(u*) 
a(u*) , 



Since S^* is bounded on L 2 (BL, C 2 ), invertible with bounded inverse, we infer from Proposition 12.51 
that a-(A LU *) = <T_(T aJ *). From (|4.15p and since the functions a and b are positive, we infer that 
there is tp £ L 2 (M) x H 2 (K), ip ^ 0, such that T^ip = Xtp for some A < and the restriction of 
Tu)*\^± > 0. From Proposition I2.6( i) conclude that U-{A U1 *) = cx_(T w ») = 1. □ 

Remark 4.5. Similarly as in Remark 14.21 we obtain, integrating traveling- wave equation (|4,lip . 

(4.17) uj* = -{F'{u)-v) a , 

where superscript a denotes average over one period, along with relation 

(4.18) g(v) a = -5u a = -5c. 

This gives the explicit co-periodic stability condition (d/dc)(F' (u*) — v*) a < 0. 
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Again, one can simplify in the homoclinic, long-wavelength limit with endstate (u, v)°° as x — > 
±00 to 

(4.19) (d/du^F'iu 00 ) - v°°) < 0, where v°° = -dg-^u 00 ). 

Comparing with the reduced traveling- wave ODE (|4.12p . we sec (cf. [PS2J) that this is the require- 
ment that (u 00 ,?; 00 ) be a saddle (center), hence by the associated existence theory we obtain again 
co-periodic instability in the homoclinic limit. In the limit towards constant states, bifurcations 
can be more complicated and both co-periodic stability and instability can occur |GFE| . 

Remark 4.6. Wolansky jW] gives a natural extension of (j4.8|) to chemotaxis systems that are "not 
in conflict" in a certain sense, with u € R m , v € R n , m and n arbitrary. For this class of systems, 
all of our results apply, giving a generalization of the stability condition to the system case. 

Remark 4.7. We mention also a more general class of chemotaxis systems identified by Horstman 
[H] . possessing a decreasing Lyapunov function £, but apparently not coercive in the sense of (jDrj) . 
It would be interesting to see whether this class of systems could be treated similarly as the case 
of viscoelasticity with strain-gradient effects, below. 

4.3. Viscoelasticity with strain-gradient effects. The equations of planar viscoelasticity with 
strain-gradient effects, written in Lagrangian coordinates, appear |BLeZ| IY1 IBYZ] as 

n - u x = 0, 

(4 20) 

u t - (VW(t)) x = {b{r)u x ) x - {h{T x )T xx ) x . 

Here r, u G R d , d = 1, 2, 3, b : R d -)■ R is a smooth function, b(r) > b > and h(a) = V 2 ^(a) > 
hold > 0, for some smooth function $ : R rf — > R. Standard choices of b and h are [BLeZl IT] 

(4.21) 6( r ) = r^ld, h{ Tx ) = Id, 

b corresponding to the usual Navier-Stokes stress tensor, and h to the phenomenological energy 
tp(p) = |p| 2 /2 of the usual Cahn-Hilliard theory of phase transitions. 

Again, we consider this system as a quasi-gradient system with constraints within the framework 
(|3.9p . where the energy is given through 



(4.22) £(t, u) := J + W(r) + (t x )) (x)dx 



and 



(4.23) M{r, u):=-[" „ ^ \ > 0. V L ,£(r, u) - ~ ^ 



K d x d x (b( T )d x )J *^"v'»^ ^ u 

The constraints of the system are given by 

C(t,u) = J r(x)dx, J u(x)dx S j . 
Further, if (r*, u*) is a solution of (|2.4p . using that the constraint function C is linear, we compute 

(4.24) iC _(-V*M)* + VH'M 0^ 

Looking for solutions of system (|2,4p . in this case we note that the n-component should be constant, 
that is u = —uj u S M. d . The r-component satisfies the equation 

(4.25) VW(t) - {VV(t x )) x + u t = 0, 

with bj r G R rf . The above equation is a nonlinear oscillator, which under appropriate conditions on 
W and \E' (for example, \P close to the identity, W not convex, or, more generally, when the matrix 

19 



pencil V 2 ^(n)/) - V 2 W(n) possesses a negative eigenvalue p) is also known to have a variety 
of solutions depending smoothly on oj t , proving that (|4.20[) has a oj = (oj t ,lo u ) E M. 2d dependent 
family of critical points, constant in the it-component. 
Note f pllBLe^l lY]) that 

d r 1 

(4.26) -S(t,u) = - {u x (x),b(r(x))u x (x))dx 

dt Jo 

is not coercive, according to the fact that \{M + M*)(t, u) = — ^ q (b(r)d )) ^ aS & ^ ar & er 
(indeed, infinite-dimensional) kernel than that of M. However, it is an interesting fact that if 
£ t £(T,u) = 0, then u x = by gjZ6]), whence d T x r t = 0, all r > 0, by (jOD]) (Tl. and so ^ 2 £{t,u) = 
\\ut\\ 2 . Thus, ^£ = gives ^(t,u) = 0, so that constant-energy solutions are stationary as in the 
coercive case (jDl'p . Likewise, we have the corresponding linear property: 

Lemma 4.8. The condition (|E1|) is satisfied at any variationally stable critical point of (|4,20p . 
Proof. Denoting Y = , assuming MA U »Y = ipY and dropping oj*, we find that 

MAY = ipY =>- Re(AY, -(M + M*)AY) = 0^d x u = 0, 

whence ipr = u x = 0, so that p = or else r = and u = constant, giving MAY = and again 
p = 0. Thus, zero is the only pure imaginary eigenvalue of MA^*. Now, suppose that Y is a 
generalized zero-eig envector, (MA) 2 Y = 0, or AM AY E kerM, i.e., d x AMAY = 0. Computing 

-(d 2 - V 2 W{r*))d x s 

-d x (d 2 x -V 2 W(r*)) d x {b{r*)d x ) 

d x AMAY = yields, in the r-coordinate, d x (d 2 — V 2 W (t*)8 x u = 0, hence, by the assumed 
definiteness of d 2 — VW(t*) on ker Tlycix*) 1 - > that d x u = 0. Substituting this into the u-coordinate 
of (|427!) then yields -d 2 (d 2 - VW(t*))t = 0, hence -d x (d 2 - VW{t*))t = 0. Combining this 
with d x u = 0, we obtain MAY = 0, and so Y is a genuine eigenfunction. This verifies that all 
zero-eigenfunctions are genuine, completing the proof. □ 

As a corollary, we obtain spectral stability whenever (|2.22p is satisfied, from which stability 
follows by the general theory developed in (YJ IBYZ1 IHoZl IKotl IZlj and references therein. Notably, 
this holds both for periodic solutions, on i7 1 [0, l] per or for front and pulse type solutions on if 1 (M). 
We note that in this case (dc/dco) has a block-diagonal structure and that (dc u /dio u ) = — l^.Thus, 
{dc/duj)(uj*) satisfies l|2321t if and only if (dc T /du T )(u}*) satisfies (f2T22l) 

In particular, following the discussion of Sturm-Liouville properties in the Calm-Halliard exam- 
ple, we conclude that for one-dimensional deformations, r, u E R , we find co-periodic stability of 
spatially periodic solutions of (|4.20[) if and only if (dc T /dui T )(uj*) > and dL/dr max > 0, that is, 
if period is increasing monotonically with amplitude. This recovers and extends the corresponding 
results of jOZlj . The phrasing of the stability condition in terms of (dc T /dca T ) yields also further 
insight, we feel, that is not evident in the original phrasing of |OZlji 10 l 



(4.27) d x AMA = d x 



For multidimensional deformations r E M. d , we obtain as a necessary and sufficient condition 
for co-periodic stability that the number of positive eigenvalues of {dc T / doj T ){uj*) is equal to the 
number of negative eigenvalues of — V 2 ^ / (r*)(9 2 + \7 2 W(t*), considered as an operator on all of 
L 2 [0, l] per (i.e., with no zero-mass restriction). This result appears to be completely new. 

Moreover, though we do not state it in detail here, we obtain a corresponding characterization 
of stability also for multiply-periodic solutions of the full d-dimensional equations of viscoelasticity 



In our notation, the necessary condition of |OZl] is effectively det(duj/du-) x det(du- /dc) > 0, with u = /(it_). 
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with strain-gradient effects, d = 1, 2, 3, for which r £ and it £ yielding d 2 + d constraints 

Cj on constrained minimization problem (|2.ip . again concerning the mechanical energy 

(4.28) £(r,it) := jf (I|u| 2 + W(t) + *(t s )) dx. 

We obtain as a necessary and sufficient condition for co-periodic stability that the number of positive 
eigenvalues of (dc T /du T )(u*), a symmetric (d 2 + d) x (d 2 + cQ-dimensional tensor, is equal to the 
number of negative eigenvalues of — V 2 ^(r*)d 2 + V 2 W(t*). For discussion of the full equations and 
their derivation, see for example [BLeZl lD]. This result, again, appears to be completely new; indeed, 
the approach by other (e.g., Evans function) means appears quite forbidding in likely complexity. 
Again, we collect the results from this subsection in the following proposition. 

Proposition 4.9. Assume that (t*,u*) is a ID-periodic or multiD periodic pattern of (|4.20p 
The following assertions hold true: 

(i) For the case of one- dimensional deformations tGR, the periodic pattern (t*,u*) is stable 
if and only if (dc T /dio T )(uj*) > and dL/dT max > ; 

(ii) For the case of multidimensional deformations r £ M. d , the periodic pattern (r*,u*) is stable 
if and only if the number of positive eigenvalues of {dc T / duj T ){uj*) is equal to the number 
of negative eigenvalues of — V 2 ^(r*)9 2 + V 2 W{t*), considered as an operator on all of 
-^ 2 [0) l]per; 

(iii) For the case of full d- dimensional equations of viscoelasticity with strain- gradient effects, the 
multiply -periodic pattern (u*,v*) is stable if and only if the number of positive eigenvalues 
of (dc T I duj T ) , a symmetric (d 2 + d) x (d 2 + d)- dimensional tensor, is equal to the number 
of negative eigenvalues of —V 2 ^{T*)d 2 + V 2 W(t*) considered as an operator on all of 

-^ 2 [0) l]per- 

Remark 4.10. It has been shown by other means in |BYZllZ3] that variational and time-evolutionary 
stability are equivalent for (14.20h : here, we recover these results as a consequence of our general 
theory. We note that it was shown in [BYZ[IZ3] also that the spectrum of MA U * is real for ReA > 0, 
which is special to the precise structure of the model, and does not hold in general. 

Remark 4.11. As in Remark 14. 2[ we find, integrating the traveling- wave ODE (|4.25p over one 
period, that uj* = — VW(T*) a , where superscript a denotes average over one period, giving the 
explicit (singly periodic) stability condition 

(4.29) d(W'(T*) a )/dc T < and dL/dr mSLX > (case d = 1) 
or a^(d(VW(T*) a )/dc T ) = a-(A w *) (case d > 1), where c T = j^T*(x)dx. 

Remark 4.12. On their natural spaces H = (H 1 x L 2 )[0, l] per and H = (L 2 x H 1 )^, l] per , respec- 
tively (the spaces on which the second variation is bounded and coercive), (|4.22p and (j4.9j) are not 
only not C 2 , but not even continuous, since t x and it enter nonlinear ly, but are not controlled in 
L°° by the norm associated with H. 

4.4. Generalized KdV. It is well-known from the celebrated work of Grillakis-Shatah-Strauss 
( |GSS| ). that the stability of solitary waves arising in models such as the Klein-Gordon equation 
or the generalized KdV equation can be studied using variational results similar to our results. We 
briefly mention here the generalized KdV model and the corresponding functionals. Consider the 
equation 

(4.30) u t + {VF(u)) x = u xxx . 

Here it £ M. m and F : W 71 — > R is a smooth function. This equation can have periodic waves 
with respect to co-periodic perturbations as shown in |BJK| or solitary waves on the whole line as 
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discussed in jGSS]. This equation also fits the general framework (|3.2p with 





In the case of periodic patterns we can take x £ [0, l] pcr and consider (|4.30p in the energy space 
H = Hp ei [0, 2], with mass constraint C(u) = (f^ u(x)dx, u 2 (x)dx). 

Remark 4.13. It is readily verified that energies (|4T2|) . P~3Tj) are C 2 on H = H 1 ^, l] per or iT^R). 
Moreover, there is a weak existence theory |Mej on the energy space H, verifying condition (|B2p . 

4.5. General conservation laws. For periodic solutions of general parabolic conservation laws 

u t + (f(u)) x = (B(u)u x ) x , 

without the assumption of generalized gradient form, Oh and Zumbrun ( [UZlj . Thm. 5.9) have ob- 
tained by direct Evans function/Stability index techniques the related necessary stability condition 

that det^(<9c/9w)(u;*)^ have a certain sign sgn7 determined by the orientation properties of u* as 

a solution of the underlying traveling- wave ODE. The Jacobian (dc/dui)(u}*) is connected to the 
(necessary) stability index condition of Theorem 5.9, |OZlj through the observation that, in the 
notation of the reference, 



where q = lo* and u* m = c. This is somewhat analogous to the relation established in |PWj between 
Evans and variational stability conditions in the Hamiltonian solitary wave case; see also |Z2j . The 
sign sgn7 may be recognized, by a homotopy argument similar to that used to prove the full index 
relation, as measuring the parity of the sign of the Morse index of the linearized traveling-wave 

operator \B(u*{x)d x + d-B(u*(x))ii*^ — df{u*{x)) equivalent (by integration) to the restriction to 
the set of zero-mass functions VC(ti*)^ of the full linearized operator 



where C(u) := j* u(x)dx, thus completing the analogy to the generalized gradient case. Indeed, 
as this argument reiterates, the stability index itself may be viewed as a mod two version of the 
classical Sturm-Liouville theory applicable to the non-selfadjoint case; see especially the discussion 
of [El [U] from a similar perspective of the general self-adjoint case, for which the zeros counted in 
the second-order scalar case are replaced by zeros of an appropriate Wronskian. 



We now describe some interesting consequences regarding the Whitham modulation equations 
and modulational stability of periodic waves, in particular, the modulational dichotomies mentioned 
in the introduction. The Whitham modulation equations, a formal WKB expansion in the high- 
frequency limit, is known (see [Sel IOZ21 iJZBl IB JNRZl IFST] and especially [NFUl INR2] ) to have a 
dispersion relation agreeing to lowest order with the spectral expansion of the critical eigenmodes of 
L u * determining low-frequency (sideband) stability. In particular, well-posedness of the Whitham 
equation may be seen to be necessary for low-frequency modulational stability. 

As critical, or slow, modes are associated to linear order with variations dX/du: and dX/ds along 
the manifold of nearby stationary solutions, the equations naturally involve the same quantities 
already studied in relation to co-periodicity, and this turns out to yield in some situations substantial 





5. Sideband instability and modulational dichotomies 
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further information. For details of the derivation of the Whitham equations in the generality 
considered here, see, for example, [Sel IOZ2j . 

5.1. Viscoelasticity with strain-gradient effects. The first-order Witham system for (|4.20p is 
(see |JNRZ2j ): 

r? ~ < = 0, 

(5.1) <- W(r)£ = 0, 

k t = 0, 

where superscript a denotes average over one period, k denotes wave number, or one over period. 
By the block-triangular structure, we may ignore the trivial third equation and concentrate on 
the first two. Indeed, this is a common feature of all of our dissipative examples, as the general 
fc-equation is k t — (<r(w, k)k) x = 0, where a denotes speed of the associated periodic traveling wave, 
which for gradient-type systems is identically zero. 
Integrating the traveling-wave ODE (12.41) 

-(V^(t x )) x + VW(t)+u = 0, 

over one period, where u is the Lagrange multiplier for the associated constrained minimization 
problem in r, we find that VW(r) a = —u, hence f)5. 1[) writes more simply as 

ct - d x = 0, 

dt + oj x = 0, 

where c = r a is the constraint for the minimization problem, d = u* is a parameter for the mini- 
mization problem (irrelevant, by Galillean invariance with respect to u), and u = co(c) independent 
of u is determined by the solution structure of the standing-wave problem. This is of exactly the 
same form as the equations of viscoelasticity themselves, whence the condition for hyperbolicity of 
the equations is that dco/dc (recall, this is a Hessian for the constrained minimization problem, so 
symmetric) be negative definite. 

It has been shown that hyperbolicity of the Whitham system is necessary for sideband stability. 
This means that periodic waves are sideband stable only ifdco/dc < 0, in which case Corollary 12.91 
implies that they are co-periodic stable only if the second variation operator 

-V 2 *(r*)^ + V 2 W{t*) 

is nonnegative on all of H l [0, l] per , not just for perturbations with zero mean. This is false in 
the scalar case r £ 1, by Sturm-Liouville theory and the fact that r* is a zero-eigenfunction that 
changes sign. Thus, we recover in straightforward fashion the dichotomy result of |OZlj . 

Theorem 5.1. In the scalar case, t,u £ R, assume that (t*,u*) is a periodic solution of (|4,20p . 
Then (r*,u*) is either co-periodic unstable or else it is sideband unstable: in all cases, modu- 
lationally unstable. For systems t,u 6 M. d , solutions are modulationally unstable unless A^* = 
-VH{T*)dl + V 2 W(t*) is nonnegative on all of H 1 [0, l]^ er . 

This yields for the first time an extension to multidimensional deformations of the dichotomy of 
[OZT] . Moreover, it gives new information even in the one-dimensional deformation case, removing 
a superfluous technical assumption (see Theorem 7.1, |OZl| ) that period T be increasing with 
respect to amplitude a. This both simplifies the analysis, eliminating the often-difficult problem 
of determining sgn(dT/da), and answers the question left open in jOZlj whether there might exist 
modulationally stable periodic waves with period decreasing with respect to amplitude. 

Further, though, as discussed in Section 14.1.11 we cannot say much in the general system case 
regarding nonnegativity of A u * , due to additional structure we can in the present situation say sub- 
stantially more. Specifically, consider the class of multidimensional deformation solutions obtained 
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by continuation from the special class of decoupled solutions, equal to a one-dimensional deforma- 
tion solution in one of its components, say, (r,-, , uj t ), and constant in the other components (rj,Uj), 
j 7^ j*. This generates a large subclass of the possible multidimensional deformation solutions: in 
some cases, perhaps, all. 

Corollary 5.2. In the (singly periodic) system case, all periodic solutions of (|4.20p lying on a 
branch of transverse solutions extending from a decoupled solution are modulationally unstable. 

Proof. By Sturm-Liouville considerations, decoupled solutions (reducing to the scalar case), satisfy 
a-(A w *) > 1, hence are modulationally unstable. Noting that Ker(A UJ *) has constant dimension so 
long as the solution remains transversal, we find that no eigenvalues can pass through the origin 
and so a-(A u] *) = 1 along the entire branch, yielding the result by Theorem [5TTJ D 

Remark 5.3. Though transversality is difficult to check analytically, it is straightforward to imple- 
ment as part of a numerical continuation study. From a practical point of view, the main importance 
of Corollary 15.21 is that waves obtained through numerical continuation from decoupled solutions 
are unstable at least up to the point that bifurcation first occurs. 

Remark 5.4. It is worth noting that the solutions bifurcating from the one-dimensional case include 
truly multi-dimensional ones and not only one-dimensional solutions in disguise. As noted implicitly 
in the description (15. ip . transverse planar periodic solutions of (|4.20j) may be locally parametrized, 
up to x-translates, by the wave number k = 1/L, where L is the period, the mean c G M d of r 
over one period, and the material velocity u = constant- the latter entering trivially, by invariance 
of (|4.20p under the Galillean transformation u — > u + c, c constant. Among these solutions are 
"disguised" one-dimensional solutions, for which r is confined to an affine set; see (Yj for further 
classification/discussion. However, a convenient asymmetry in the coefficient 1/T3 of the parabolic 
term in u for the standard choice of viscosity (|4.2ip forces such solutions to lie either in the 
hyperplane T3 = 1 or else along the line t\ = ■ ■ ■ = r^-i = 0, so that solutions can be effectively 
one-dimensional in this sense only if T3 = 1, or else t\ = ti = 0. But, these possibilities can be 
excluded by choosing a r-mean c for which C3 7^ 1 and (ci, . . . , C4-1) 7^ (0, . . . , 0), of which there 
are uncountably many. 

Remark 5.5. Theorem 15.11 resolves a question left open in [OZlj whether there might exist modu- 
lationally stable waves in the case dL/dr mSiX < 0. The reason for the difference in these two sets of 
results is that the treatment of co-periodic stability in [OZlj proceeds by a stability index (Evans 
function) computation, which gives information on the parity of the number of unstable eigenvalues 
of the linearized operator about the wave and not the actual number as determined here. 

5.1.1. The multiply periodic case. It is interesting to consider the implications of sideband stability 
for multiply periodic solutions of the full d-dimensional equations of viscoelasticity with strain- 
gradient effects, for which the first-order Whitham homogenized system is (ignoring the decoupled 
/c-equation) 



x G M d , with c now a matrix G W dxd and d G R d , or cu — V x (V x u;(c)) = 0. For this system, 
well-posedness (hyperbolicity) is well-known [D] to be equivalent to rank-one convexity of (dc/du), 
a less-restrictive condition than the convexity condition arising in the one-dimensional case x G M 1 . 
(Recall, under this notation, c = r a G M rfxa! .) This enforces only d positive eigenvalues in the 
signature of (dc/du>), with another d(d— 1) zero eigenvalues forced by rotational symmetry, leaving 
d 2 — d possible negative eigenvalues, in principle compatible with variational stability. At the 
same time, the relation between variational and time-periodic stability is unclear, as the argument 



(5.3) 



c t - V x d 
dt + V x u 
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of Lemma 13.191 yielding (|Eip breaks down in the fully multidimensional case0 Thus, it seems 
possible that one might find stable multiply periodic solutions of the equations of viscoelasticity. 

5.2. Cahn Hilliard systems. The second-order Whitham modulated system for (|4.ip may be 
seen by a similar (straightforward) derivation as in [Sej to be 

, . u a t =-d 2 M(u)-VW{u))\ 

1 ' ' h = d x {a{u\ k)d x (u a ) + P(u a , k)8 x k). 

where a, f3, since irrelevant for our considerations here, are left unspecified. Integrating the 
traveling-wave ODE (12.4j) 

-u" + VW(u) +u = 0, 

over one period, where oj is the Lagrange multiplier for the associated constrained minimization 
problem, we find that (VW(u) — u") a = X7W(u) a = —oj, hence (|5.4p writes more simply as 

^ ^ c t = -8%(u) = -d x ((doj/dc)d x c), 

kt = d x (a(c, k)d x c + /3(c, k)d x k). 

where c = u a is the constraint for the minimization problem. This system is well-posed (parabolic), 
corresponding to sideband stability, only if doj/dc is negative definite, with the same conclusions 
as in the case of viscoelasticity just considered. 

In particular, in the scalar case, recovering a result of |Holj . we may immediately conclude 
existence of a dichotomy as in jOZlj . stating that waves are either sideband unstable or co-periodic 
unstable; in either case, unstable. As in the previous subsection, we obtain also extensions to the 
system case. 

It is worth mentioning that the one-dimensional theory just described extends essentially un- 
changed to the multiply periodic case, i.e., stationary solutions in 12" [0, l]^ er . Again, there is 
a single constraint c equal to the mean over [0, l] m , and a single Lagrange multiplier oj, with 
standing-wave PDE 

Au - VW(u) = -oj. 

In the scalar case, the condition for co-periodic stability is dc/doj > 0, by our abstract theory. 
On the other hand, the Whitham modulation equations are easily found to be 

ct = Aw = — V x ■ ((doj/ dc)V x c), 

h = V x - (a(c, k)V x c + /3(c, k)V x k), 

k £ R m , hence are well-posed only if (doj/dc) < 0. We have therefore that multiply periodic waves 
are either co-periodically unstable or else sideband unstable, a dichotomy similar to that observed 
for singly periodic scalar waves, implying that multiply-periodic stationary solutions of the scalar 
Cahn-Hilliard equation are always modulationally unstable under assumption (|H3p , So far as we 
know, this result for multiply-periodic Cahn-Hilliard waves is new (though see Remark 15.71 below) ; 
the connection to sideband stability likewise appears to be a novel addition to the Cahn-Hilliard 
literature. 

In the system case, we obtain a partial analog, similarly as for viscoelasticity, namely, co- 
periodically stable solutions of (|4.ip are necessarily sideband unstable, unless A^* = — A x + 
V 2 W(u*) is nonnegative on H^O, l]™ er - We collect the results of this subsection in the follow- 
ing theorem. 

Theorem 5.6. Assume u* is a periodic solution of (|4.ip . Then u* is either co-periodic unstable 
or sideband unstable. In the system case, co-periodically stable waves are sideband unstable unless 
A w * = -A x + V 2 W(u*) is nonnegative on L 2 [0, l]™ r . 



Specifically, [M + M*) is no longer definite in the u coordinate, a key element in the proof. 
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Remark 5.7. In the scalar, one-dimensional case, we may obtain the result of modulational (but not 
necessarily sideband) instability more directly, as in |GNj . by the observation that, by Sturm- 
Liouville considerations, periodic waves cannot be co-periodically stable on a doubled domain 
[0,2T], since the unconstrained problem has at least two unstable modes, only one of which may 
be stabilized by the constraint of constant mass. Recalling that modulational stability is equiva- 
lent to co-periodic stability on multiple periods (0, NT) for arbitrary N, we obtain the result. In 
the two-dimensional case, we may obtain, similarly, co-periodic instability on the multiple domain 
(0, NT] x (0, NT], by the Nodal Domain Theorem of Courant |CHj . using a bit of planar topological 
reasoning to conclude that there must either be a nodal domain in each cell that is disconnected 
from the boundary, or else a nodal curve which disconnects two opposing edges of the cell bound- 
ary, and thus there must be at least 2N nodal domains in (0, NT] x (0, NT], and thereby 2N — 1 
unstable modes, again yielding modulational instability for N > 2. This argument is suggestive 
also in dimensions m = 3 and higher, but would require further work to eliminate the possibility 
that there are only two nodal domains in (0, NT] m , for example a thickened m-dimensional lattice 
(no longer disconnecting space in dimensions m > 3) and its complement. 

Remark 5.8. It would be interesting to consider the consequences of our theory in the planar multi- 
dimensional case considered by Howard in [Ho2j . that is, the implications for transverse instability. 

5.3. Coupled conservative-reaction systems. For equations (|4.8p . it is readily seen that the 
second-order Whitham modulated system is 

u a t =d x (b(u) a d x (F'(u)-v) 
(5.6) ) 

k t = dja(u a , k)8 x (u a ) + p(u a , k)8 x k ) , 



so that, denoting c = u a , b(u) a = (3(c,k) > 0, and referring to the traveling-wave equation (|4.1ip . 
we obtain the equation 

ct = -((3(c,k)uj x ) x = -(/3(c,k)(duj/dc)c x ) x , 
k t = (q(c, k)c x + f3(c, k)k x ) x , 

which is well-posed (parabolic) only if (dc/duj)(uj*) < 0, in contradiction with co-periodic stability. 
Thus, we find a dichotomy as in the previous cases. 
Likewise, for multiply periodic solutions of 

^ 5 ^ ut = Vjr • (a(u)V x u - b(u) V x vj , 

v t = A x v + 5u + g(v), 
or stationary solutions in -ff 1 [0, l]™ er , we obtain 

c t = -Vje • ((3(c,k)d(u/dc)V x c), 

kt = y x - (a(c, k)V x c + P(c, k)V x k), 

k G M m , where c is equal to the mean over [0, l] m , and oj is the corresponding Lagrange multiplier, 
again yielding a modulational dichotomy unless A^, is not nonnegative on i/ 1 [0, l]™ er , We collect 
the results of this subsection in the following theorem. 

Theorem 5.9. Assume that (u*,v*) is a periodic solution of (|4.8[) . Then (u*,v*) is either co- 
periodic unstable or else it is sideband unstable: in all cases, modulationally unstable. Multiply 
periodic solutions of (15 ,7p are modulationally unstable unless A^* is nonnegative on H 1 ^, 1]™. 

5.4. Generalized KdV. Results of [BJ21j for the generalized KdV equation, relating sideband 
stability to Jacobian determinants of action variables, appear to be an instance of a similar phe- 
nomenon in the Hamiltonian case. This would be a an interesting direction for further study. 
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5.5. Dichotomies revisited. The above examples can be understood in a simpler and more uni- 
fied way through the underlying quasi-gradient structure of the linearized generator C = —MA. 
Namely, we have only to recall that sideband instability is defined as stability for f 6l, |£| << 1 
of the small eigenvalues of the Floquet operator 

(5.8) £ € := e-^ x Ce^ x = M^A $ 

bifurcating from zero eigenvalues of Cq = C, where the £ subscript indicates conjugation by the 
multiplication operator e l ^ x , or 

Ct. := e-^Ce^, := e'^Me^, A^ := e -* x Ae* x . 

Modulational stability, similarly, is defined as stability for £ G R, of M^A^. 

Note first that A^ inherits automatically the self-adjoint property of A, and also the quasi- 
gradient property 

(5.9) i(M € + M*) > for all |£| « 1, 

as seen by the computation (M 5 + M|) = (M + M*) € = (M + M*))J 2 (M + M*) 1 ^ 2 , both of which 
follow from the adjoint rule 

(5.10) (Ntf = (JV*) e> 

where the * on the lefthand side refers to L 2 [0, L] pcT and the * on the righthand side refers to L 2 (R) 
adjoint, valid for any operator N that has an L 2 (R) adjoint N*, is L-periodic in the sense that 

(5.11) (Nf)(x + L) = (Nf(- + L))(x) 

for C°° test functions /, and is bounded from W s '°° to L°° for s sufficiently large: in particular, for 
the periodic-coefficient differential operators considered here in the applications. We note further 
that, under these assumptions, and N£ take .L-periodic functions to L-periodic functions. 

The latter statement follows by noting that N% is also invariant under shifts by one period. The 
former follows by noting that, for u, v periodic, H s and satisfying 

(5.12) u{x + L) = 7ti(x), v(x + L) = jv(x), \ j\ = 1, 

we have that (u, Nv)xj2\Qn is equal to L times the mean of u ■ Nv over R. This follows by smooth 
truncation on bigger and bigger domains, with truncation occurring over a single period, and noting 
that this bounded error goes to zero in the limit. But, then, this is also equal to the mean of N*u-v 
over R, by taking the L 2 adjoint for the truncated approximants and arguing as before by continuity. 
Finally, observe that (e~ l ^ x Ne l ^ x u, v) = (Nu, v) , where u = e l ^ x u and v = e l ^ x v satisfy the above 
property ()5.12p . Thus, it is equal to (u,N*v) = (u, (N*)^v) as claimed. For periodic-coefficient 
differential operators p(d x ) = Sj=o a j^-> (|5.10j) . and the property that periodic functions are taken 
to periodic functions, may be verified directly, using the rule 

(p(&c))$ =P{dx + ii)- 

Now, assume, as is easily verified in each of the cases considered (but need not always be true), 
that 

(Fl) kerM^ = 0, 

so that the corresponding linear evolution equation (d/dt)Y = —MAY is unconstrained (has no 
associated conservation laws). Then, for M self-adjoint, we obtain the correspondence 

(5.13) a -{A) < a^(A^) = n_(M^),EB 



The strict inequality cr-(A) < cr-(Af) is possible in the cases we consider, due to bifurcation of ker A. 
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where n_ is the number of negative eigenvalues counted by algebraic multiplicity, as may be seen by 

— 1/ 2 1/2 1/2 1/2 

the similarity transformation -)• M $ ' M^M*' = M£ 1 A^M^ (Here, we are assuming 

sufficient regularity in coefficients that eigenfunctions remain in proper spaces upon application of 

±1/2 

and also implicitly that spectrum is discrete, properties again easily verifiable in each of 
the cases previously considered). More generally, the result carries over to the case (|Df '|) provided 
(M e + M s ) > 0; see Remark [3J2D 

Thus, we find immediately that modulational stability is violated unless is unconditionally 
stable for |£| << f , which, in the limit £ — > 0, gives unconditional (neutral) stability of A. Moreover, 
again invoking continuity of spectra with respect to £, we find that, assuming co-periodic stability, 
or conditional stability of A, that the negative eigenvalues of M^Ae must be small eigenvalues 
bifurcating from zero-eigenvalues of MA, which, by the property (|Eip as guaranteed by Lemma 
13. 19^ are the only neutral (i.e., zero real part) eigenvalues of MA. 

This argument shows that, assuming the property (|5.13p . first, modulational stability requires 
unconstrained stability of A, and, second, if co-periodic stability holds, then sideband stability 
requires unconstrained stability of A. The second observation may be recognized as exactly the 
modulational dichotomy just proved case-by-case, the first as a generalization of the one obtained 
by nodal domain consideration in Remark 15.71 

Remark 5.10. We obtain by this argument modulational dichotomies for Cahn-Hilliard and coupled 
conservative-reaction diffusion equations. On the other hand, we don't (quite) obtain a result 
for the viscoelastic case, for which M is not self-adjoint, and so our simple argument for (|5. 131) 
does not apply. We can recover this however by a "vanishing viscosity" argument, noting for 
e > that n_(M(e)^4 f ) = 0-(A £ ) for M(e) = M + eld and, by property (jPTjl . n (M(e)^) = 
no(M(s)^A^) = ao(A^) = 0. Mimicking the proof of Lemma 14.81 we find, similarly, that the center 
subspace of M^A^ consists of the zero-subspace of M^A^, which is likewise empty. Thus, no zeros 
cross (or reach) the imaginary axis during the homotopy and we obtain the result in the limit as 
e — > 0, for any £ ^ sufficiently small that A^ does not have a kernel. 

Remark 5.11. The above clarifies somewhat the mechanism behind observed modulational di- 
chotomy results; namely, the origins of such dichotomies are a priori knowledge of modulational 
instability (encoded, through considerations as above, or as in Remark 15.71 in unconstrained insta- 
bility of A Wt ) plus the central property (|Eip that the center subspace of MA Wf consists entirely of 
ker MA^ . 

6. Discussion 

We conclude by briefly discussing how our results connect to other stability studies. 

6.1. Relation to the Evans function. First, we note the following connection to the Melnikov 
integral and the Evans function: 

Remark 6.1. If r = 1, the quantity 

(dc/du) = (Vc, (dx/doj)) = -{A^dx/d^^dx/duj)) 

is a Melnikov integral involving variations in uj along the manifold of nearby stationary solutions, 
and as pointed out in |PWllZ2] and elsewhere should correspond to the first nonvanishing derivative 
of an associated Evans function. 

Of course, we don't need this connection, since we have an if and only if condition for asymptotic 
stability under our assumption (|Eip . which is perhaps an analogy in the dissipative context to the 
[GSSj assumption in the Hamiltonian setting that J be one to one (the difference being that f|El|) 
seems to be satisfied for all the systems we know of/are interested in, while J one-to-one does not 
hold for (gKdV) and other primary examples). 
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6.2. Localized structures: problems on the whole line. The results of |HoH IPSH IZ3j on 

exponential instability of solitary, or "pulse-type" exponentially localized spikes solutions suggest 
another dichotomy different from the one we proved in Section [5j if essential spectrum of the 
linearization L = —MA W * along the spike is good, then point spectrum must be bad. The two 
dichotomies are apparently different, but do resemble each other — for, co-periodic stability is 
related to point spectrum, sideband to essential spectrum (thinking of the homoclinic limit [G])F^I 

Moreover, we can recover (slightly weakened versions of) the above mentioned results on insta- 
bility of spikes from [HoU IPS11 IZ3| using our results on variational stability of periodic solutions. 
Below we will illustrate the main ideas in the specific case of the viscoelasticity model (|4.20p with 
d = 1. First, we assume that the essential spectrum of the linearization along the spike is stable, 
otherwise the instability is trivial. We note that the quasi-gradient structure implies that it is 
enough to prove variational instability in order to obtain (neutral) time-evolutionary instability. In 
\Z3\ lY] it was shown that for the viscoelasticity with strain-gradient effects model (|4.2U|) the two 
types of instability are actually equivalent. 

We suspect this is true for a larger class of models, but choose not to pursue this here. Rather, 
we just notice that instability of A^* , without constraint implies already an exponentially localized 
initial data for the evolution problem on the whole line for which £ is nondecreasing in time and 
initially strictly less that £(X*), whereas the manifold of nearby solitary wave solutions, consisting 
entirely of translates of X*, has energy = £{X*). Thus, by continuity, X(t) cannot converge in H 2 
to the set of translates of X*, else the energies would also converge; this shows that the standard 
notion (see |Z3[ IT]) of L 1 n H s — > H s orbital asymptotic stability cannot hold for s > 2. 

Next, we remark that any homoclinic spike solution can be approximated by a sequence of 
periodic patterns having large period 2T. This is true for any system that has a Hamiltonian 
structure, in particular the models ()4. 1|) . (|4.8p and (|4.20p . Then, we use the convergence results for 
the point spectrum from [Gj ISS| to prove the (above type of neutral) instability of the homoclinic 
spike solution by showing that the approximating periodic pattern is unstable for any T > T* > 0. 

Averaging in (|4.25p . one readily checks that 

(6.1) w T = --L f W'(T*(x))dx -> -W'(t°°), as T -»• oo. 

2T J_ T 

Similarly, using the definition of the constraint function, we have that 

(6.2) ^- = -^-[ T*(x)dx^T QO , as T^oo. 
K J 2T 2T J_ T y ' 

From (I6.ip and (I6.2p we conclude that 

(6.3) T*^- ->■ —W"(t°°), as T ->■ oo0 

OC T 

Since r* is a solution of (|4.25p having finite limit at ±oo we infer that W"(t°°) > 0. Our claim 
follows immediately from Proposition 14.91 (i) and (|6.3|) . yielding a sequence of eigenvalues ut > 
converging as T — > oo to an eigenvalue v of A^* . With further work, exploiting the results in 
|SSj . one may show, provided the limiting operator A^* has a spectral gap at v = 0, that the 
associated eigenfunctions Jt have uniform exponential decay at ±oo, hence extract a subsequence 
converging to an eigenfunction / of A u * ; moreover, the same uniform exponential decay allows us to 
conclude that / is both orthogonal to the translational eigenfunction X* and satisfies the linearized 
constraint, with eigenvalue v > corresponding to "neutral" (i.e., nonstrict) orbital instability as 

In view of the discussion of Section 15.51 the link may be just that modulation and the large-period limit arc 
related through small-£, or "long-wave" response, both having the effect of "removing" constraints. 

14 As the reverese limit dc/dui — » oo indicates, solutions with different multipliers to are infinitely far apart in the 
whole-line problem, and do not play a role; in particular, ker MA = Span{X*} is one-dimensional. 
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claimed. Provided that Ker(A UJ *) is spanned by the translational mode X*, we obtain in fact strict 
variational instability v > 0. 

On the other hand, we have, more simply, just by Sturm-Liouville considerations, that A is 
variationally unstable without constraints, without the above construction. Either argument leads 
to the conclusion that asymptotic orbital stability cannot hold. If one can verify further that the 
derivative of the standard Evans function does not vanish at A = 0, then one could go further 
(by property (El)) to conclude strict, or exponential instability as shown in [Z3| (and analogously 
for reaction diffusion-conservation law spikes in |PS2j ). Such a result would then yield by [G] 
the known result of co-periodic instability for sufficiently large period, completing the circle of 
arguments. 

Remark 6.2. As noted in [Z3] , stability of front-type solutions may also be studied variationally for 
viscoelasticity with strain-gradient effects. It would be interesting to try to phrase this entirely in 
terms of the structure ([3.90 : in particular, we suspect that the similar result of [LWj in the context 
of chemotaxis is another face of the same basic mechanism. 

Remark 6.3. In the cases considered in [HolllP"ST}IZ3| . the existence problem is scalar second order 
with spectral gap at v = 0, whence it is easily by dimensionality seen that X* is the only element (up 
to constant multiple) of Ker( J 4 aj *), hence we indeed recover strict constrained variational instability 
from the limiting periodic argument above. However, we do not see a correspondingly easy way 
to see a priori that the derivative of the Evans function does not vanish at the origin, so for the 
moment obtain by this argument only neutral instability for the time-evolutionary problem. 

Remark 6.4. The above instability arguments (including those of [Holt IPST1 IZ3j ) rely on spectral 
gap of Aj* at the origin, v = 0, and the associated property of exponential convergence of X* to its 
endstate as x — > dboo. In the critical case of an algebraically-decaying homoclinic solution, it has 
been shown for the Keller-Segel equation that stability holds, with algebraic decay rate [CFIIBCC] . 

6.3. Coupled conservative reaction diffusion spikes on bounded domains. As pointed 
out in |PS2j . coupled conservation-reaction diffusion spikes have been numerically observed \Dr\ 
IHHMOl IKSl IKR] . suggesting that they may sometimes be stable on finite domains. As we pointed 
out earlier, the stability with periodic boundary conditions is equivalent to stability with Neumann 
boundary conditions, since odd eigenfunctions are automatically stable, due to the fact that the 
translational derivative possesses only two sign changes in a domain of minimal period. Assuming 
monotonicity of the period, one concludes that patterns are stable whenever dc/duj > 0. Now letting 
the period tend to infinity, there are just a few possibilities for a periodic solution: convergence to 
a homoclinic solution, convergence to a heteroclinic solution, or unbounded amplitude. In fact, all 
scenarios occur in relevant circumstances. We already discussed the case of the spike limit, in which 
all periodic solutions are unstable. In the Cahn-Hilliard example with, say, a cubic nonlinearity 
W'(u) = —u + ii 3 , one finds stable solutions periodic solutions that converge to a heteroclinic loop 
as the period tends to infinity. In fact, the energy is coercive in this case, so that there necessarily 
exists a constrained minimizer for any prescribed mass. For masses \C(u)\ < l/y/3, homo geneous 
equilibria are unstable so that coercivity of the energy enforces the existence of stable equilibria 
for all periods, which, given a priori bounds from the energy, therefore necessarily converge to a 
heteroclinic loop (or a pair of layers /kinks) . Similar considerations apply to the coupled conservative 
reaction-diffusion context when we have a globally coercive energy. In the case of chemotaxis, this 
global coercivity fails. Also, the nonlinear oscillator that describes periodic solutions possesses only 
two equilibria, which excludes heteroclinic loops. In the large- wavelength limit, one finds spike-like 
solutions which do not converge to a homoclinic solutions, but diverge along a family of spikes to 
infinity as the period goes to infinity, consistent with our previous discussion. 
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